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Abstract 

Possible effects of large-scale magnetic fields on the Cosmic Microwave Background 
(CMB) are reviewed. Depending on the specific branch of the spectrum of plasma ex- 
citations, magnetic fields are treated either within a two-fluid plasma description or 
within an effective (one-fluid) approach. The uniform field approximation is contrasted 
with the fully inhomogcncous field approximation. It is argued that the interplay be- 
tween CMB physics and large-scale magnetic fields will represent a rather interesting 
cross-disciplinary arena along the next few years. 
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1 Why CMB anisotropies could be magnetized? 



Simplified magneto-hydrodynamical estimates imply that the magnetic diffusivity length 
scale^ in the interstellar medium is of the order of the astronomical unit. On the other 
hand, magnetic fields are present over much larger length-scales so that there seems to 
be rather compelling evidence that galaxies, clusters and possibly super-clusters are all 
magnetized. The question that arises naturally in this context concerns the possible effects 
of large-scale magnetic fields on the Cosmic Microwave Background (CMB). In the last fifty 
years various cosmological mechanisms for the origin of large scale magnetic fields have been 
proposed. While different mechanisms rely on diverse physical assumptions some general 
features can be identified: 

• the majority of the cosmological mechanisms imply the existence of large-scale mag- 
netic fields after equality (but before decoupling); 

• the magnetic field present after equality is, according to the mentioned mechanisms, 
fully inhomogeneous; 

• the typical amplitudes and length-scales of the magnetic field are characteristic of the 
given model. 

The first possibility we can think of implies that magnetic fields are produced, at a 

given epoch in the life of the Universe, inside the Hubble radius, for instance by a phase 

transition or by any other phenomenon able to generate a charge separation and, ultimately, 

an electric current. In this context, the correlation scale of the field is much smaller that 

the typical scale of the gravitational collapse of the proto-galaxy which is of the order of 

the Mpc. In fact, if the Universe is decelerating and if the correlation scale evolves as the 

scale factor, the Hubble radius grows much faster than the correlation scale. Of course, one 

might invoke the possibility that the correlation scale of the magnetic field evolves more 

rapidly than the scale factor. A well founded physical rationale for this occurrence is what 

is normally called inverse cascade, i.e. the possibility that magnetic (as well as kinetic) 

energy density is transferred from small to large scales. This implies, in real space, that 

(highly energetic) small scale magnetic domains may coalesce to form magnetic domains of 

smaller energy but over larger scales. In the best of all possible situations, i.e. when inverse 

cascade is very effective, it seems rather hard to justify a growth of the correlation scale 

that would eventually end up into a Mpc scale at the onset of gravitational collapse. In 

Fig. ^ we report a schematic illustration of the evolution of the Hubble radius Rh and of 

^The magnetic diffusivity length is the typical scale below which magnetic fields are diffused because of 
the finite value of the conductivity of the (interstellar) medium. 
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INSIDE THE HUBBLE RADIUS 




Figure 1: Evolution of the correlation scale for magnetic fields produced inside the Hubble 
radius. The horizontal thick dashed line marks the end of the radiation-dominated phase 
and the onset of the matter-dominated phase. The horizontal thin dashed line marks the 
moment of e~^~e~ annihilation (see also footnoote 2). The full (vertical) lines represent the 
evolution of the Hubble radius during the different stages of the life of the Universe. The 
dashed (vertical) lines illustrate the evolution of the correlation scale of the magnetic fields. 
In the absence of inverse cascade the evolution of the correlation scale is given by the (inner) 
vertical dashed lines. If inverse cascade takes place the evolution of the correlation scale is 
faster than the first power of the scale factor (for instance a^/^) but always slower than the 
Hubble radius. 
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the correlation scale of the magnetic field as a function of the scale factor. In Fig. ^the 
horizontal dashed line simply marks the end of the radiation-dominated phase and the onset 
of the matter dominated phase: while above the dashed line the Hubble radius evolves as 
(where a is the scale factor), below the dashed line the Hubble radius evolves as a^/^. 

We consider, for simplicity, a magnetic field whose typical correlation scale is as large 
as the Hubble radius at the electro-weak epoch when the temperature of the plasma was 
of the order of 100 GeV. If the correlation scale evolves as the scale factor, the Hubble 
radius at the electroweak epoch (roughly 3 cm) projects today over a scale of the order 
of the astronomical unit. If inverse cascades are invoked, the correlation scale may grow, 
depending on the specific features of the cascade, up to 100 A.U. or even up to 100 pc. In 
both cases the final scale is too small if compared with the typical scale of the gravitational 
collapse of the proto-galaxy. In Fig. ^a particular model for the evolution of the correlation 
scale A(o) has been reported ^. 

In the context of (conventional or unconventional) inflationary models, on the contrary, 
the correlation scale of the produced magnetic fields may be rather high. The physical 
picture is that during the inflationary phase the quantum fluctuations of the hypercharge 
field are amplified either through some coupling to the geometry or through the coupling to 
some other spectator field. Since during inflation the Hubble radius is roughly constant (see 
Fig. [21), the correlation scale evolves much faster than the Hubble radius itself and, therefore, 
large scale magnetic domains can naturally be obtained. Notice that, in Fig. |21the (vertical) 
dashed lines illustrate the evolution of the Hubble radius (that is roughly constant during 
inflation) while the full line denotes the evolution of the correlation scale. Furthermore, 
the horizontal (dashed) lines mark, from top to bottom, the end of the inflationary phase 
and the onset of the matter-dominated phase. This phenomenon can be understood as 
the gauge counterpart of the super-adiabatic amplification of the scalar and tensor modes 
of the geometry. The main problem, in such a framework, is to get large amplitudes for 
scale of the order of the Mpc at the onset of gravitational collapse. Models where the 
gauge couplings are effectively dynamical (breaking, consequently, the Weyl invariance of 
the evolution equations of Abelian gauge modes) may provide rather intense magnetic fields. 

The two extreme possibilities mentioned above may be sometimes combined. For in- 
stance, it can happen that magnetic fields are produced by super-adiabatic amplification 
of vacuum fluctuations during an inflationary stage of expansion. After exiting the hori- 
zon, the gauge modes will reenter at different moments all along the radiation and matter 

^Notice, as it will be discussed later, that the inverse cascade lasts, in principle, only down to the time of 

e+ e~ annihilation (see also thin dashed horizontal line in Fig. since for temperatures smaller than 

Tg+_g- the Reynolds number drops below 1. This is the result of the sudden drop in the number of charged 
particles that leads to a rather long mean free path for the photons. 
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dominated epochs. The spectrum of the primordial gauge fields after reentry will not only 
be determined by the amplification mechanism but also on the plasma effects. As soon 
as the magnetic inhomogeneities reenter, some other physical process, taking place inside 
the Hubble radius, may be triggered by the presence of large scale magnetic fields. An 
example, in this context, is the production of topologically non-trivial configurations of the 
hypercharge field (hypermagnetic knots) from a stochastic background of hypercharge fields 
with vanishing helicity. 

Since magnetic fields are produced over different length scales, it seems plausible to 
discuss their effect on the anisotropics of the Cosmic Microwave Background radiation 
(CMB). The finite values of the thermal and magnetic diffusivity scales forbids the presence 
of magnetic fields over length-scales that arc very small in comparison with the size of the 
Hubble radius at a given epoch. On the other hand, dissipation is not effective in erasing 
the magnetic energy density over very large scales comparable with the Hubble radius itself. 

Different topics are relevant for the interplay between large-scale magnetic fields and 
CMB physics. A simplified list of the main effects involving, simultaneously, large-scale 
magnetic fields and CMB physics must include: 

• possible distortions of the CMB blackbody spectrum; 

• possible induction of primary anisotropics; 

• possible effects on the CMB polarization. 

To this list it is appropriate to add another important topic that has to do with the sub- 
traction of foregrounds. In fact, to disentangle the signal of the CMB anisotropics one has 
usually to subtract the emission of the galaxy. The synchrotron emission is indeed one of 
the foregrounds to be subtracted. 

To estimate the effects of large-scale magnetic fields on CMB physics various approx- 
imations have been used through the years. Again, for sake of simplicity, it is useful to 
divide them into two broad classes: 

• uniform field approximation; 

• inhomogeneous field approximation. 

In the uniform field approximation, the magnetic field is taken to be uniform and pointing 
along a specific direction in the sky. In this approximation (pioneered by Zeldovich in the 
sixties) the magnetic field slightly breaks the spatial isotropy of the background geometry. 
The break in the spatial isotropy can be easily connected with a difference in the propagation 
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SUPERADIABATIC AMPLIFICATION 




Itf'cm 



Figure 2: Evolution of the correlation scale if magnetic fields would be produced by supera- 
diabatic amplification during a conventional inflationary phase. The dashed vertical lines 
denote, in the present figure, the evolution of the Hubble radius while the full line denotes 
the evolution of the correlation scale (typically selected to smaller than the Hubble radius 
during inflation). 

of electromagnetic signals along different spatial directions. This phenomenon is ultimately 
responsible, in this context, for the anisotropy of the CMB. 

In a complementary perspective, fully inhomogeneous large-scale magnetic fields do not 
break the spatial isotropy of the background geometry. This happens, for instance, when 
the magnetic field is fully inhomogeneous and characterized by its two-point function. In the 
uniform field approximation the resulting evolution equations are rather tractable and semi- 
analytical estimates of various effects can be swiftly obtained. It is, however, problematic 
to imagine that a uniform magnetic field (oriented along a specific spatial direction) may 
be generated in the early Universe. The presence of a uniform field could be viewed as a 
remnant of the initial conditions of the evolution of the Universe. 

Conversely, the fully inhomogenous approach is more realistic but also harder to im- 
plement analytically. An intermediate approach is to discuss the uniform field case sup- 
plemented by the appropriate fluctuations. In this approach the parameters describing the 
model are given by the uniform component of the magnetic field and by the two-point 
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function of the fluctuations. 

An important point to be borne in mind when discussing the effects of large-scale mag- 
netic fields on CMB anisotropies concerns the suitable set of equations used to model the 
plasma effects. For the purposes of this paper we will confine ourselves to the two most 
practical sets of equations normally employed in this type of studies, namely 

• a set of two-fluid plasma equations in curved space; 

• an effective one-fluid (magnetohydrodynamical) description. 

These two systems of equations are not equivalent for all the frequencies of the spectrum 
of plasma excitations. In particular, if we ought to address phenomena involving the prop- 
agation of electromagnetic disturbances in the plasma (as in the case of Faraday) rotation, 
the two fluid description (or even a kinetic Vlasov-Landau description) have to be em- 
ployed. Magnetohydrodynamics (MHD) is a good effective description for sufficiently small 
frequencies in the spectrum of plasma excitations. 

The plan of the present review is the following. In Section 2 the uniform field approx- 
imation will be reviewed. In Section 3 a two-fluid plasma description will be introduced. 
As an application, the Faraday effect in the uniform field approximation will be discussed. 
In Section 4, after discussing a one-fluid plasma description, fully inhomogeneous magnetic 
fields will be analyzed. Section 5 contains an introduction to the evolution of metric fluctu- 
ations in the presence of large-scale magnetic fields. In this context, the problem of initial 
conditions of CMB anisotropies will be reviewed in the presence of large-scale magnetic 
fields. The final part of Section 5 contains a summary of the numerical results obtained so 
far in the context of fully inhomogeneous magnetic fields. A brief collection of concluding 
remarks and future perspectives is collected in Section 6. 

2 Uniform field approximation 

The uniform field approximation has a long history that could be traced back to the work 
of Zeldovich (see also |HE1)- Before the formulation of inflationary models the problem 
of the the initial conditions of Friedmann-Robertson- Walker models was debated. Among 
different options (early dominance of viscous effects, anisotropy in the expansion, etc.) 
a particularly interesting possibility was represented by a uniform background magnetic 
field oriented along a specific spatial direction. Indeed, close to a big-bang singularity the 
spatial gradients are sub- leading but the model can well by anisotropic @1 El • In other 
words, close to the big bang, the homogeneous (but anisotropic) metric can fall in one of 
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Figure 3: Uniform field approximation. A simplified geometrical configuration where the 
magnetic field is directed along the x axis. 

the Bianchi classes [7? and a uniform magnetic field can fit into this scheme. During an 
inflationary phase the anisotropy decays exponentially and, therefore, we have to imagine 
that the (uniform) field considered in this section was created, for instance, close to the end 
of inflation. Since the work of Zeldovich PQE]) the uniform field approximation has been 
exploited in various frameworks (see Refs. |H1 El El El E] and references therein for a 
tentative list of theoretical studies). 

Consider then, as an illustration, the case when the magnetic field is oriented along the 
X axis. Suppose also, for sake of simplicity, that the Universe is filled by a perfect barotropic 
fluid and that the spatial curvature vanishes. In this case the line element must account for 
the possibility of a different expansion rate along the x axis or along the (y, z) plane and 
may have the form 

ds'^ = dt^ - a^(t)dx^ - h^(t){dy^ + dz'^). (2.1) 

The metric 1)2. 1() corresponds to a particular case of a Bianchi-type I metric. Uniform 
magnetic fields can be introduced also in other Bianchi metrics P^ll() | [TT l ll2j. The anisotropy 
in the expansion and the total expansion rate are, in this case, 

. H-F ^ H + 2F , , 

A = — e = ^^, (2.2) 

where H = d/a and F = b/b are the expansion rates and the overdot denotes a derivation 
with respect to the cosmic time coordinate t. 

In the model defined by Eq. ()2.1|1 the propagation of an electromagnetic signal will be 
different along the x and y axes (see Fig. Therefore, the temperature anisotropy will 
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depend both on G and A: 

^ = ^^^-lAedt, (2.3) 

where we used that = To/a and Ty = T^/h together with the identities = e~ -f 
and 6-1 = 6-/^'='*. 

The quantities appearing in Eq. (|2.3p can be estimated by studying the evolution of 
the anisotropic expansion and of the total expansion rate along the different stages of the 
life of the Universe. For this purpose it is practical to phrase Einstein equations directly 
in terms of A and Q. The procedure is rather simple since it suffices to write the extrinsic 
curvature in terms of H and F. Then, by recalling the definitions of Eq. 1)2. 2(1 . the full set 
of Einstein equations and of the continuity equations can be written, respectively, as 

"2-3e 



+ 



(2.4) 



2 

Ae + A{e + 3e^) = -i6TTGpq, (2.5) 

p + e(4 + 3e)p = 0, (2.6) 

q = 3qee+^q@A, (2.7) 



where the following conventions have been adopted 



w = e+-, q= — . (2.8 
3 p 

In Eq. 1)2. 8jl w is the barotropic index of the fiuid background and 

PB = — , /9B+4F/)B=0. (2.9) 
ovr 

In the case e = the barotropic index is exactly the one of radiation and Eqs. ()2.4() - ()2.7() 
admit a solution where ^ ~ and @ ~ {2t)~^. Inserting then Eq. 1)2. 5(1 into Eq. 1)2. 7|) (and 
recahing that 6 = -29^) it is easy to obtain the evolution of q{t), i.e. 

lit) = / X , t > h. (2.10) 



l + 4giln (^^ 

This means that in the case of a radiation background the anisotropy in the expansion is 
conserved and the ratio between radiation and magnetic energy densities is roughly constant 
(up to logarithmic corrections). For generalization of this analysis to more complicated 
situations involving different equations of state or various fluids and fields see, for instance, 

[Hiiiniiini- 

The property expressed by Eq. (|2.10() is a consequence of the fact that both energy- 
momentum tensors (i.e. the one of the background fluid and the one of the magnetic field) 
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are traceless. Therefore, the evolution of q{t) will be different if the background fluid is 
dominated by dusty sources. In the case of dusty matter Eq. 1)2. 8(1 implies that e = —1/3. 
Therefore, an approximate solution of the system in this case is given by 

m^Y^, ^(t)=.-12g,q(^), q(t)=q,^^^y (2.11) 

By following the evolution of the system through equality up to decoupling the value of the 
anisotropy can then be roughly calculated to be 

AT ^2 



Zdec (2.12) 



By requiring that the temperature anisotropy is smaller than 10 ^ we get in the case of 
^^tot = 1, 

Bo<2x IQ-^ Gauss, (2.13) 

where Bq refers to the value of the field at the present time. The logic leading to the limit 
reported in Eq. ()2.13|) has been followed, for instance, in Ref. 

A complementary way of describing the effect of a uniform magnetic field on the temper- 
ature anisotropies relies on the observation that, thanks to the breaking of spatial isotropy, 
the angular power spectrum will not depend only on the polar angle but also on the az- 
imuthal angle. In fact, if spatial isotropy is unbroken, the angular power spectrum is 
customarily defined as 

C{^) = (Ai(ni,ro)Ai(n2,ro)), (2.14) 

where Aj is the brightness perturbation in the intensity of the radiation field whose expan- 
sion in terms of spherical harmonics is 

oo e 

Ai(n,ro) = ^ aumYimin). (2.15) 

1=0 m=-i 

If the background space is isotropic, the ensemble average of a£m only depends upon i, not 
upon m, i.e. 

{aemae'm') = CeSu'Smm'- (2-16) 

Therefore, inserting Eq. (|2.15|) into Eq. (|2.14jl . the addition theorem of spherical harmonics, 
i.e. 

Piini ■ h2) = yi*mini)YUn2), (2.17) 

m=—£ 

leads to the usual form of the angular power spectrum: 

C(7?) = (Ai(ni, ro)Ai(n2, ro)) ^ ^ ^(2£ + 1)QP,(t?). (2.18) 
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where 'd = fii • n2- 

If spatial isotropy is slightly broken (as in the case when a uniform magnetic field is 
present), Eq. (|2.16j) is not justified since the ensemble average of the aim should also depend 
upon m and not only upon i. Indeed, it can be shown that when a uniform magnetic field 
is present, the average of temperature fluctuations in two directions of the sky can be 
expressed as 

(Ai(ni,ro)A(n2,ro)) = ^ Q(m)Y£,m(ni)y^*„(n2) 

+ Y.Di{m)[Ye+i,^{hi)Y;_,^^{h2)+Ye.i,^{M^^^ (2.19) 

where both Ce{m) and Di{m) depend upon the (uniform) magnetic field intensity. From 
Eq. (|2.19|) it is apparent that the presence of a magnetic field with uniform component 
induces off-diagonal correlations in multipole space between a given £ and multipoles £±2. 
The breaking of spatial isotropy in the CMB maps can be constrained, for instance, from 
WMAP data . 

This analysis has been recently performed in jl4j (see also ^H]). The constraints on 
the off-diagonal correlations in multipole space imply a bound on the (present) uniform 
component of the magnetic field intensity 

Bo < 1.7 X 10"^ G. (2.20) 

It should be mentioned that the bound of Eq. (|2.2()j) has been derived by taking into account 
also the fluctuations of the magnetic field around the uniform configuration. The presence 
of the fluctuations of a uniform configuration introduces extra parameters (like the spectral 
index of the fluctuations) that make the allowed region in the parameter space effectively 
two-dimensional. Various studies assume a (uniform) background field and consider the 
possible fluctuations of such a configuration (see for instance ^1 El Ej and references 
therein). In particular, in Ref. the magnetic field has been parametrized in terms of 
a uniform component directed along a specific axis and a tangled component. The Alfen- 
waves modes associated with this configuration may induce a small rotational perturbation 
in the last scattering surface. 

In this analysis has been performed in a non-relativistic approximation where the 
scalar perturbations of the geometry obey linearized Newtonian equations of motion. The 
magnetic field has been assumed to be uniform (but supplemented with the appropriate 
fiuctuations). In this framework, it has been argued that the effect of the presence of the 
(uniform) magnetic can be phrased into an effective renormalization of the speed of sound 
of the baryons (see also J20*). 
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3 Two-fluid plasma description and CMB polarization 



Right before decoupling the value of the Debye screening length is much larger than the 
electron and photon mean free paths. In fact, the Debye length scale can be written as 



Ad = ~ 10 f .""Q ^ ( -^^] cm. (3.1) 

where no — n-g ~ n\ is the mean electron-ion density and Tgi is the electron-ion temperature. 
For length scales larger than Ad, the plasma is globally neutral. For typical values of the 
parameters in Eq. (|3.1|) . Ad is of the order of 10 cm around recombination while the electron 
and photon mean free paths are much larger, i.e. 

4 5.7 X lO^cm, ~ 10''(1 + z)-2(17b/io)"^ Mpc, (3.2) 

where z = ao/a — 1 is the redshift, fib the baryon fraction of the critical energy density and 
/iQ the indetermination in the Hubble parameter. 

The other set of quantities relevant for the determination of the appropriate system of 
equations is given by the plasma and Larmor frequencies of electrons and ions, i.e. 



Up 



" ^ ( 10^ ) '"^^ " ^^-^H iF^ ' ^^-^^ 



rip. 



Note that the typical value of B around decoupling reported in Eqs. (|3.3|) and (|3.4jl corre- 
sponds, today, to a magnetic field of the order of 10~^-10~^^ G. 

If we are now interested in the propagation of electromagnetic disturbances with fre- 
quencies larger than the plasma and Larmor frequencies, then it is clear that a two-fluid 
plasma description should be adopted. For some classes of problems it will also be useful 
to employ directly a kinetic (Vlasov-Landau [^1211123]) description. If, on the contrary, 
we would be interested in the spectrum of plasma excitations for w < Wp, then effective 
(one- fluid) descriptions (like MHD j24j l can be usefully employed. This second regime will 
be addressed in Section 4. 

The two-fluid plasma equations in curved space-time will now be discussed and treated 
in analogy with the well-known flat-space case Let us consider, for simplicity a con- 
formally flat geometry of Friedmann-Robertson- Walker (FRW) type: 

ds^ = a^{T)[dT'^ -dx\ (3.5) 

The following two equations account for the continuity of electrons and ions charge densities 

ng + SweTCue + (we + 1)V • (ngiJg) = 0, (3-6) 

n[ + SwiTirii + {wi + 1)V • (riiVi) = 0, (3.7) 



11 



where the prime denotes a derivation with respect to the conformal time coordinate r while 
the other quantities appearing in Eqs. H3.6() and 1)3 .7(1 are 

TL = a' /a^ rie = a'^fie, ni = a^hi. (3-8) 

In Eqs. (|3.6|) and 1)3. 7|) and Wi are the barotropic indices for the electron and ion fluids. 
Both electrons and ions are non-relativistic (for the ion mass, in the numerical estimates, 
we take m[ ~ nip where nip denotes the proton mass) . Hence the barotropic indices We and 
Wi will be close to zero. In fact, the energy and pressure densities of an ideal electronic gas 
are given by 

Pe = nc(me + -Tc], Pe = neTe, (3.9) 



2 

and since We,! = Te,i/me,i, we will also have We,! <C 1 as far as Te,i <C m^ie- Recalling that 
Pe = rUe n-e and pi = mini , the momentum conservation equations are 



Pe[Ve + + {Vc ' ^H] = -rieci E+—xB], (3.10) 



Pi[v'i+nvi + ivi-V)vi] = nie\^E + ^ X Bj, (3.11) 

where E = a^£ and B = a^B are the conformally rescaled electromagnetic fields obeying 
the following set of equations: 

V-^ = 47re(ni-ne), (3.12) 

V-B = 0, V X E = --B', (3.13) 

c 

^ ^ 1 ^ 47re 

V X B = -E' H {niVi - n^Ve). (3.14) 

c c 

The dispersion relations relevant for the propagation of electromagnetic waves with uo > uop 
can be found, after some algebra, in the case of a uniform magnetic field directed, for 
instance, along the z axis^. 

Equations 1)3. 6(1 and (|3.7|1 can be first linearized around a globally neutral configuration, 

i.e. 

Tie (r, x) = no + 5ne(r, x), ni(r, x) = no + (5ni(r, x), (3.15) 

where no is the common electron and ion density. The remaining equations can be also 
linearized around the uniformly magnetized configuration: 

B{r],x) = Bo + 6B{'n,x), 

Ve,i{r],x) = 6ve,ii'n,x), E{t],x) = 5E{rj,x). (3.16) 



*In principle we should allow a different expansion rate along the (x, y) plane and along the z axis, as 
discussed in the previous Section (see, for instance, Eq. 12.11 1. This is not strictly necessary for the derivation 
of the dispersion relations. While the anisotropic expansion caused by a uniform magnetic field affects the 
primary anisotropics, the dispersion relations are weakly affected by the curved space corrections. This 
statement will be explicitly verified in our case. 
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The linearized equations around the uniformly magnetized background can be solved j26j . 
The dielectric tensor can then be obtained and it has the following general form j26j : 



where 



e{uj,a) = —ie2{u!,a) ei(u;,a) , (3-17) 




^2 ^2 

6||(^,a) = l- - (3.18) 

^gi(« + l) ^ge(« + l) 

ei w, a = 1 2? — r"n:2 2 2/ , ^^2 2-' (3.19) 

u;^(a + 1)^ - a;^(a + 1)^ - (jgg 

£2(0;, a) = 2? — r~n2 — 2 27 — r~r\2 — 2-- (^-^o) 



The quantity ey (w, a) denotes the dielectric constant parallel to the direction of the uniform 
component of the magnetic field intensity while ei(w,a) and €2(00, a) denote the dielectric 
constants determined by the motion of charged particles orthogonal to the direction of the 
uniform magnetic field. If the uniform magnetic field vanishes, we can obtain the dielectric 
tensor of the cold plasma by formally taking the limit €2(00,0) ^ and ei{uj, a) — > ey {to, a). 
It is then clear that in the absence of magnetic field the only scales entering the dispersion 
relations are the plasma frequencies of electrons and ions. The quantity a appearing in 
the dielectric tensor is defined as a = iTC/u and is the lowest order correction arising in 
curved space. This correction is negligible for the present purposes since it is of the order 
of the ratio between the Hubble radius at the decoupling and the inverse of the frequency 
of propagation. Recalling that H^cc — 5.6 x 10^2(f2p/i2-j-i/2 ^^i^t u » 100 MHz, we 

have that \a\ < 10~22. 

To derive the dispersion relations it is appropriate to fix the coordinate system: for 
instance we can set = and ky = ksinO, = kcosO with the magnetic field oriented 
along the z direction. Introducing the refraction index n and recalling that k = oj/v = nu /c 
(where v is the phase velocity) the dispersion relations are determined from the solution of 
the following equation: 

1 1/1 1 



\eii 71-^ 



v? 2 Ul ^ CR. 



where c{9) = cos 6 and s{9) = sinO. In Eq. 1)3.21(1 the following linear combinations have 
been also introduced: 

^2 ^2 

eR = ei+62 = l- . , J' 7- r . J." 1, (3.22) 

, ,2 2 

1 Pe pi /o r,o\ 

= ei - £2 = i r—, — — T 7 r—, — — T— y [6.26) 
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The quantity eR and el are the dielectric constants appropriate for the propagation of the 
two (left and right) circular polarizations. The dispersion relations for the propagation of 
electromagnetic disturbances parallel (i.e. 6 = 0) and orthogonal (i.e. 6 = tt/2) to the 
uniform component of the magnetic field can then be obtained from Eq. 1)3. 21() : 

(n2-eR)(n2-eL) =0, ^ = 0, (3.24) 

(n2-e||)[n2(eL + eR)-2eLeR] =0, 9 = ^. (3.25) 

Equation (|3.24|) gives the usual dispersion relations for the two circular polarizations of 
the electromagnetic wave, i.e. = eR and = eL, while Eq. ()3.25|) gives those for the 
"ordinary" (i.e. n? = e||) and "extraordinary" (i.e. = 2eReL/(eR + eL)) plasma waves. 

Consider now an electromagnetic wave that is linearly polarized. A linearly polarized 
electromagnetic wave can be always seen as the sum of two circularly polarized waves. Now 
the left and right polarizations have different phase velocities and, therefore, after the waves 
have travelled a distance Az in a (cold) magnetized plasma the resulting phase difference 
will be 

Az. (3.26) 
Equation ()3.26)) has been derived using two physical assumptions 

• the plasma is cold, i.e. the temperature of the electrons does not enter the dispersion 
relations; 

• there are no other sources leading to a rotation of linearly polarized radiation. 

Let us now see what happens by relaxing these two assumptions. The adoption of a cold 
(as opposed to warm) plasma description means that the temperature of the electrons 
has been taken to be effectively zero. The effects arising from the finite temperature of 
the electrons do not modify the leading result obtained in the context of the cold plasma 
theory. The derivation of the dispersion relations in the case of a warm plasma can be 
performed, for instance, within a kinetic approach where it can be shown, following the 
same calculation discussed in the flat space case :27_ , that the first correction to the leading 
cold plasma calculation can be recast in an effective redefinition of the plasma frequency 
for the electrons, namely ujpe — > where 7 = (1 — (f^))~^/^ and (v^) is the thermal 

average of the electron velocity. 

Equation ()3.26|) has been derived under the assumption that the only source of birefrin- 
gence is represented by a uniform magnetic field. This is certainly legitimate in a standard 
electromagnetic context. There could be however other sources of cosmological birefrin- 
gence leading to a generalized expression of Faraday rotation. In particular let us consider 
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the case when, in the plasma, a dynamical pseudo-scalar field is present together with a 
magnetic field. This situation may easily arise in the case of a quintessential background 
with pseudo-scalar coupling to electromagnetism described by an action 



\g^^d^adpa - W{a) + -^^F^pF^P 



(3.27) 



where F^p and F"^ are the Maxwell field strength and its dual. The coupling of electro- 
magnetism to the field a changes the form of Eqs. 1)3. 12() . ()3.13|1 and ()3.14|1 . The resulting 
system of equations can be discussed in the presence both of a uniform magnetic field and 
of a cold (or warm) plasma. This calculation has been performed in Ref. j^H]- While Eqs. 
1)3. 13|) do not change, the generalization of Eqs. (|3.12() and (|3.14|) becomes: 



V X i? 



47re(ni - n^) + ^Va ■ B, 



-E' 



£ 

M 



a'B + VaxE 



H (niVi 

c 



UeVe). 



(3.28) 
(3.29) 



To understand the rationale of the physical difference induced by the presence of a consider, 
indeed, the case when Vcj = and the charge and current densities are vanishing. In this 
case, taking the curl of Eq. (|3.29j) and using the second of Eqs. (|3.13j) we obtain the 
following equation 

B" - c^V'^B + £ca'V X B = 0. (3.30) 



M 



Equation 1)3. 3U() implies that the two polarizations of the magnetic field are mixed as a 
consequence of the time-dependence of a. If the plasma contribution is included, Eq. (|3.21|) 
is modified '26^ . Consequently, the dispersion relations given in Eqs. 1)3. 24() and (|3.25|) are 



also modified: 
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(3.31) 
(3.32) 



The presence of the quintessence field introduces a further frequency scale into the problem, 
namely uj^ = c{(3 / M)a' . If the quintessence field dominates today (or between redshifts 
and 3), cr(to) ~ h? /m where, typically, A ~ 10"'^ eV, and m ~ o" can 

be estimated, and it is, today, (T(to) — rn^ / Hq. Hence, recalling that prior to quintessential 
dominance, d scales as a~^, from the previous expressions 



Ha 



/ M 



X 10"*^ Hz 



(3.33) 



where the values of A and m are the ones discussed above with M ~ Mp. The value of /3 
is rather uncertain and a conservative limit from radio-astronomical analyses would imply 
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/?^10"3 for 0^ z ^1 |2HlIini- Notice that the action (|in7|l can also be relevant in a class 
of baryogenesis models jSUl |^ (see also IHl] ) ■ 

From Eq. (|,3.3H) the generalized Faraday rotation experienced by the linearly polarized 
CMB travelling parallel to the magnetic field direction can be obtained as 

This expression includes, simultaneously, the birefringence produced by the evolution of 
the pseudo-scalar field a and the rotation induced by the finite value of the magnetic field 
intensity. In Ref. j34j the effect of parity-violating interactions on CMB polarization has 
been recently discussed in the absence of a magnetic field. A similar analysis was presented 
in [Hlj (see also (HE])- In the absence of magnetic field, Eq. (|3.34jl leads to a A$ that is 
frequency independent. On the contrary, if the contribution of the pseudo-scalar field a is 
negligible, then A$ decreases at high frequencies as l/w^. In the intermediate situation 
more complicated dependences of A$ upon the frequency may be envisaged. This aspect 
answers the second question raised after Eq. H3.26() : on top of the magnetic fields there 
are different sources affecting A<1>. However, the characteristic frequency dependence of 
the different signals may allow, at least in principle, to disentangle different sources of 
birefringence. 



3.1 Magnetized Faraday effect at decoupling 

Let us now assume, for simplicity, that the only source for A<^ is represented by a uniform 
magnetic field. Then, recalling the hierarchy between the plasma, Larmor and propagation 
frequencies, i.e. Eqs. (|3.3j) and (|3.4|1 . Eq. (|3.26|1 can be expanded for 



\^p,c,i/oj\ <. 1, |a;B,e,i/'^| < 1- 



(3.35) 



Thus, from Eq. 1)3. 26(1 . we will have: 
A$ 
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Az, 



(3.36) 



where, by virtue of Eqs. 1)3. 3|) and 1)3. 4|1 the sub-leading contribution of the ions has been ne- 
glected. In some simplified treatments this approximation is made from the very beginning 
by defining two generalized refraction indices tir and n\^: 

,2 



p,e 



u;(a;±a;B, e) 



(3.37) 
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( the plus and minus signs correspond, respectively, to R and L^) 

A last expansion for |u;b,c/'^| ^ 1 and for lajp^c/"^! ^ 1 brings Eq. (|3.36|) to the standard 
form, i.e. 



2c 



^ ^^^u^^^^Az. (3.38) 



From Eq. (|3.38() . by choosing a generic propagation direction, we can also derive the Faraday 
rotation rate, i.e. 

where the inonization fraction Xe has been introduced by replacing uq XeriQ and where 
we used, from Eq. (|3.38|) . that dz = cdt. In Eq. (|3.39|) . B • h denotes the projection 
of the (uniform) magnetic field intensity along the direction of propagation of the linearly 
polarized radiation with frequency u = uj/2it; t is, as usual, the cosmic time coordinate 
related to the conformal time coordinate r as a{T)dT = dt. 

For some applications it is sometimes useful to write Eq. ()3.39|) in apparently different 
forms. Recalling the definition of differential optical depth 

e = XeUoaT— = — — . (3.40) 

ao z + l 

Defining with tq the conformal time at which the signal is received, the optical depth will 
then be 

/■TO Q 

e(T, To)=/ XenoaT—dr (3.41) 
Jt ao 

while the visibility function ^ is 

/C(r) = e'e-^(^'^o). (3.42) 

Equation 1)3.42(1 gives the probability that a CMB photon was last scattered between r and 
T + dr (usually in the literature /C(r) is denoted with gir)). 

Using Eq. ()3.41|) . Eq. (|3.39() can be written as 

where the explicit expression of the Thompson cross section in terms of the classical radius 
of the electron has been used, i.e. 

2 



8vr 9 e / . ,x 

^T = —rl ro = 2. 3.44 

3 mpC"' 



^Recalling that ui = fec/riR, l, by compute du/dk we can easily get, from Eq. 113.3711 . and approximate 
expression of the group velocity. 

® Recently the possible modifications induced by large-scale magnetic fields on the thermal history of the 
Universe during the "dark age" beteween z ~ 1000 and z ~ 10 has been studied in |33| . 
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Equation (|3.43|) can be also written in terms of the typical values of the magnetic field 
intensity and of the frequency of emission. Denoting with v the physical frequency, i.e. v = 
z^/a, it is clear that BjiP' = B jv^ since the quantities in the numerator and in denominator 
redshift in the same way. Thus Eq. 1)3. 43() becomes: 

S), = _ = 3.56 X 10- (^) (— — ) (3.45) 

The typical quantities reported in Eq. 1)3. 45() refer to the 'present time. It is also conventional 
to use the same quantities but evaluated at the decoupling epoch. In this case the magnetic 
field will be of the order of the mG while the physical frequency (at the decoupling time) 
will be of the order of 100 THz. 



3.2 Faraday effect and radiative transfer equations 

The Faraday rotation rate derived in Eq. (|3.43|) directly enters the radiative transfer equa- 
tions. From Eq. (|3.39jl it is possible to give an order of magnitude estimate of <I> itself. 
Indeed, by integrating Eq. ()3.39|) with respect to i, we can obtain that 

e" ^ f 

$ = 7^ B - h Xe nodt. (3.46) 

2 TT c i/^ J 

Recalling that a{T)dT = dt^ the integral appearing in Eq. I|3.46j) can be easily estimated 
using Eq. (|3.41() : the result will be 

Xe nodt = lillM. (3.47) 

For orders of magnitude estimates, one can safely assume that the optical depth, i.e. e(r, tq) 
is of order unity out to redshift of decoupling, when the polarization is generated j3 711381 OU] . 
Therefore, after averaging over the polatizations of the magnetic field, the final result can 
be expressed in terms of \/ {^'^): 

B \/30GHz\2 



To obtain more accurate estimates of the Faraday rotation rate, the relevant radiative 
transfer equations must be integrated [HZl (see also HDI for an introduction to the radiative 
transfer equations in CMB physics). As it is well known, the polarization of the CMB is 
a higher order effect with respect to the temperature anisotropies and it depends crucially 
upon the properties of the photon phase space distribution that should have a non- vanishing 
quadrupole moment. Indeed, to zeroth order in the tight-coupling expansion photons and 
baryons are synchronized so well that the CMB is not polarized. This means, technically. 
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that the brightness perturbation associated with the Q Stokes parameter (i.e. Aq) is zero 
to zeroth order in the tight couphng expansion. Since Aq measures the degree of linear 
polarization, the CMB, to this order, is not polarized. However, to first order in the tight- 
coupling expansion a small amount of linear polarization is generated and it is proportional 
to the expansion parameter, i.e. 1/e', and to the (zeroth-order) dipole. 

This argument disregards the important fact that the Q and U Stokes parameters (and 
their associated brightness fluctuations Aq and Au) are not invariant under rotations. It 
is possible to take into account this aspect by defining a generalized degree of polarization 
that is proportional to (A^ + Aq). 

The presence of a magnetic field introduces then further terms in the radiative transfer 
equations. These terms couple directly the evolution of Aq to the evolution of A\j and the 
strength of the coupling is given by the Faraday rotation rate 0;$ introduced in Eq. (|3.45|) . 
Therefore, the relevant radiative transfer equations can be written as ^ 

A'q + {ikfi + e')AQ = |[1 - P2ip)]SQ + 2J7$Au, (3.49) 

A'y + {ikn + e')Au = -2n^AQ, (3.50) 

where, 

Sq = A12 + Aqo + Aq2. (3.51) 

Eqs. (|3.49|) and (|3.5fl|) are derived by using the following expansions for the brightness 
perturbations: 

AQ(fc, h, r) = Y^i-^y + 1) ^Q^(^' ^) Pii^ ■ (3-52) 

Av{k, n, t) = Y.i-iY (2^ + 1) ^veik, r) Pe{k ■ n), (3.53) 

e 

where k is the momentum of the Fourier expansion, k its direction; h is the direction of the 
photon momentum and fi = k ■ h. Consequently, in Eq. ()3.51|) . A12 denotes the quadrupole 
moment of the brightness fluctuation associated with the I Stokes parameter, while Aqq and 
Aq2 denote, respectively, the monopole and quadrupole moment of Aq. For completeness, 
and to anticipate further considerations on the cross power spectra, we also introduce the 
expansion for the brightness perturbation associated with the intensity of the radiation field: 

Ai(^, n, r) = ^(-i)^ {2£ + 1) Aj^fc, r) P^(fe • n). (3.54) 

e 

An important remark concerns here the symmetry of the radiative transfer equations. 

In Eqs. H3.49() and 1)3. 5U() /i denotes the projection of the Fourier mode along the direction 

^Here we are only discussing the Faraday effect induced on the polarization in the case of scalar fluctua- 
tions of the geometry. 
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of the photon momentum, i.e. fi = k ■ h. In Eqs. p. 49(1 and 1)3. 5U() the magnetic field 
may introduce an extra preferred direction, so that the evolution equations will depend, 
ultimately, not only upon /i but also upon B -n (see Eq. (|3.43jl ). Furthermore, the analysis 
may become even more cumbersome if the magnetic field is allowed not to be uniform (as 
assumed in the present and in the previous section). To make Eqs. ()3.49|) and (|3.5U|) more 
tractable the direction of the magnetic field intensity is taken to lie along the direction of 
the Fourier mode. Instead of working with Aq and Au, it is useful to write down Eqs. 
1)3. 49(1 and 1)3. 5U() directly in terms of the two associated complex conjugate combinations, 
i.e. 

7W± = AQ±iAu. (3.55) 
Using Eq. (|3.55|) . Eqs. (|3.49|) and 1)3.50(1 give, after some simple algebra, 

M'± + {ikfi + e'± 2in^)M± = ^e'(l - h^)Sq. (3.56) 

Using the definitions given in Eqs. p.41|l and ()3.42|) . Eq. I|3.56() can be integrated easily 
and the result of the integration is 

4 Jo 

where is given by Eq. (|3.45|) . This equation can be used for semi-analytical estimates of 
the mixing between Aq and Ay to corroborate numerical estimates of the same effect. 

The analysis performed in terms of Aq and Au can also be expressed in terms of the 
E and B modes. In the latter language, the mixing between Aq and Au translates into a 
mixing between E and B modes ^. 

Recalling, in fact, that under clockwise rotations of if A4± transform as 

M± = e^'^''PM± (3.58) 

the combinations 

M±{n) = ^a±2/m ±2yim{n) (3.59) 
can be expanded in terms of the spin-2 spherical harmonics, i.e. ±23^™ ('^)5 with 

a±2/m = Jdn ±23^£m(^) (^Q ± «^u)(^). (3.60) 

The "electric" and "magnetic" components of polarization are eigenstates of parity and may 
be defined as 

^em = ~2 + «-2/m)i = -(a2/m — 0--2/m)- (3.61) 

*It should be clear that the E abd B modes discussed here are not necessarily related to physical electric 
and magnetic fields. 
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These newly defined variables are expanded in terms of ordinary spherical harmonics, 

Yem{n), 

E(n) = J2 afmYUn), B(n) = J] af^Yemin). (3.62) 

£m Im 

In connection with the spin-2 spherical harmonics appearing in Eqs. (|3.59|) and (|3.60() it 
is relevant to mention here that a generic spin-s spherical harmonic is a known concept in 
the quantum mechanical theory of angular momentum (see, for instance, |41j). A typical 
quantum mechanical problem is to look for the representations of the operator specifying 
three-dimensional rotations, i.e. R] this problem is usually approached within the so-called 
Wigner matrix elements, i.e. l^l^^/iR) = (ji m > where j denotes the eigenvalue 

of and m denotes the eigenvalue of J^. Now, if we replace m' — > — s, j £, we have the 
definition of spin-s spherical harmonics in terms of the "D^l ^{a, /3, 0), i.e. 



^/^^p2,^(a,/?,0), (3.63) 



s3^^m(a,/3) 

where a, (5 and 7 (set to zero in the above definition) are the Euler angles defined as in 
[IT] . If s = 0, 2?^^^^(a,/3,0) = v/(2^ + l)/47rY£m(a,/5) where Yim{ot,l3) are the ordinary 
spherical harmonics. 

To define properly the cross power spectra we also recall the expansion of the intensity 
fluctuations of the radiation field (defined in Eq. I|3.54() ). in terms of spherical harmonics: 

Ai(n)=^aLy,^(n), (3.64) 

Im 

where we wrote explicitly since the fluctuations in the intensity of the radiation field, 
i.e. Ai are nothing but the fluctuations in the CMB temperature. 

Under parity inversion, the components appearing in Eqs. 1)3. 62j) and (|3.64|1 transform 

as 

O^m ~^ a^mi «Fm ~^ ^^mi ^fm ~^ ^J^. (3.65) 

Therefore, the E-modes have the same parity of the T-modes which have, in turn, the same 
parity of spherical harmonics, i.e. (—1)^. On the contrary, the B-modes have (— l)^"^^ parity. 

The existence of linear polarization allows for 6 different cross power spectra to be 
determined, in principle, from data that measure the full temperature and polarization 
anisotropy information. The cross power spectra can be defined in terms of the spectral 
functions ' where X and Y stand for E, B or T depending on the cross-correlation one 
is interested in: 

Cf^y k^dk y (3.66) 
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Therefore, if we are interested in the TE correlations we just have to set X = T and 
y = E and use the relevant expansions given above. In the following, we will denote the 
correlations as TT, EE, BB, TB and so on. This notation refers to the definition given in 

Eq. (Hmn) . 

Let us now see how the definition p.66() works. Suppose we are interested in the TT 
correlations, i.e. the usual and well known temperature correlations. From Eq. ()3.66|) we 
will have 

Cr = ^Jk^dk± ^""-g^;f-^'^ (3.67) 

m=—l ^ ' 

Now, from Eq. (|3.64|) . using the orthogonality of spherical harmonics, we have that 

aJmif^) = d'>T'Yem{n)Ai{k,h). (3.68) 



Inserting Eq. ()3.68|) into Eq. 1)3. 67|) and recalling the expansion of Ai(A;,n) in terms of 
Legendre polynomials we get 

qTT f dkk'^\Au\'^. (3.69) 
vr J 

To get to Eq. p.69p the following two identities have been used, i.e. 



/ 



47r 



- 2^ + 1 - 

J2 Yl^{k)YUn) = -^Piik ■ n). (3.70) 

m=-t 

The second identity in Eq. 1)3. 70(1 has been already exploited in Eq. (|2.17jl . 

In similar ways, different expressions for the other correlations may be obtained. Notice 
that Eq. 1)3. 69() is a consequence of the specific conventions adopted in Eq. (|3.54() . In 
particular note that the a factor [21 + 1) is included in the expansion. It must be clearly 
said that this is matter of conventions. For instance, in Refs. [IH] and [121 the factor (2£ + l) 
as well as the factor {—iY are included in the expansion ^. On the contrary, in Ref. the 
authors did not include the factor (2^ + 1) in the expansions for the brightness perturbations 
that we defined in Eqs. (|3.52|) . (|3.53j) and (|3.54|) . Consequently, the cross correlations will 
inherit extra-terms that are simply a consequence of the different conventions adopted. So, 
for instance, in Eq. (|3.69|) . and in the conventions of [Si; a factor (2^ + 1)^ typically appears 
in the denominator. 



^Notice that the authors of Ref. )39| do not include the factor (— i)* but do include the factor (2^+ 1) in 
the expansion of the various brightness perturbations. 
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Recalling the connection between the Wigner matrix elements and the spin-s spherical 
harmonics, i.e. Eq. ()3.63|1 . it is possible to show that, under complex conjugation, 

«i2/n^ = (-l)"^«T2/n^, (ajf'^)* = (3-71) 

where the second equality follows from the first one by using Eq. (|3.61() . It is then possible 
to show that while the TB and EB correlators are parity-odd, all the other correlators (i.e. 
TT, BB, EE, TE) are parity even. 

3.3 Limits on uniform magnetic fields from Faraday effect 

In the absence of magnetic fields, the scalar fluctuations of the geometry only generate 
E-modes. The presence of a magnetic field induces, however, a mixing between E-modes 
and B-modes. Consequently, if we have an initial E-mode generated by scalar metric per- 
turbations, the presence of a magnetic field induces a non-vanishing B-mode. In the case 
of scalar fluctuations, a non-vanishing TE power spectrum arises directly as a result of the 
CMB polarization as it can be verified by studying the first-order tight coupling expansion. 
Recently, the WMAP collaboration [441 145j measured the TE correlations. As discussed in 
j46j . the induced TB angular power spectrum is, roughly 0^ times smaller than the initial 
TE correlations. From the WMAP measurement, the measured value of the TE correlation 
and the absence of the TB correlations allows to set a limit on il^. At a given frequency, 
the limit on can be translated into a bound on the magnetic field intensity, i.e. 

Bo < G (3.72) 

at a typical (present) frequency of 30 GHz ^6j. The values of the TE and EE correlations 
induced by the presence of a magnetic field are $7^ times smaller than their values in the 
absence of the magnetic fields. In Fig. ^ the autocorretaion of the i?-modes of CMB 
polarization is reported in the case of a uniform magnetic field (recall that 0^ is defined in 
Eq. ()3.45|) '). Notice that when 217^ ~ 1 the TE (as well as the the other correlations) are 
strongly damped and the differential Faraday rotation has a net depolarizing effect j39j . 

We conclude this section by remarking that Faraday rotation is not the only way to 
generate i?-modes. In fact: 

• B-modes can be generated directly by fully inhomogeneous magnetic fields (see fol- 
lowing sections); 

• B-modes can be generated by gravitational waves (in Fig. Q the GW contribution is 
reported for a value of the tensor-to scalar ratio compatible with the present WMAP 
data); 



23 



e-17 



10 



100 



1000 



I 



Figure 4: Angular autocorrelation of B-modes of CMB polarization generated by a uniform 
magnetic field in a cosmic concordance model (adapted from Ref. 46 ). Notice that the 
vertical axis of has been rescaled by {2^^)~'^ (see Eq. (|3.45|) '). With the dashed line marked 
with "lensing" the possible B-mode arising from weak lensing is reported. The B-mode 
induced by GW is also reported (dotted line) for a typical value of the tensor to scalar ratio 
r = 0.01. 
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• B-modes can be generated by weak lensing. 

To this list one should also add the possibility that some other form of birefringence is 
allowed (as elaborated in the derivation of Eq. (|3.34|) '). In connection with the GW contri- 
bution, it should be noted that r (denoting the ration between the tensor and scalar (pri- 
mordial) power spectra) is limited from above, by current experimental data, i.e. r < 0.53. 
In Fig. |l]the contribution of the B-mode induced by GW is illustrated in the case r = 0.01. 

The peculiar feature of the Faraday rotation signal induced by a uniform magnetic field 
is the frequency dependence. In fact, by going at higher frequencies, the signal gets reduced 
as u~^. On the contrary the B-modes induced by gravitational waves does not depend upon 
the frequency. The birefringence induced by a dynamical quintessence (pseudo-scalar) field 
is also independent on the frequency. The need of multifrequency determination of Faraday 
rotation is then clear and this point has been also correctly emphasized in [471 148j . In |17j a 
scheme for the treatment of Faraday rotation in the case of fully inhomogeneous magnetic 
fields has been studied and along a similar line another analysis has also been presented in 

In the following years on the CMB side the Planck mission 50 will give us, hopefully, 
important clues concerning CMB polarization. For the purposes of the present subject it is 
appropriate to recall the Planck will operate in 9 frequency channels ranging from 30 to 857 
GHz. Multifrequency Faraday rotation measurements may become possible. Unfortunately, 
owing to the fact that the nominal Planck sensitivity to polarization is larger at higher 
frequencies, these measurements may also be difficult. In fact, Faraday rotation signals are 
suppressed by the square of the frequency and are therefore smaller at high frequencies. 
Other instruments dedicated to the analysis of CMB polarization may also give important 
clues. These experiments include polarization-sensitive balloon borne experiments (like the 
polarization-sensitive version of Boomerang j51[ I52j and its possible improvements) as well 
as ground based arrays like the promising QUIET (Q/UImagingExperimenT) j53j . On the 
radio-astronomical side, the very ambitious SKA (Square Kilometer Array) project will be 
able to give us precious full-sky surveys of Faraday rotation measurements (31] (see also 
j55[ I56j l. The SKA experiment should cover a frequency range between 0.1 and 20 GHz. 

4 One-fluid plasma description 

Since a fully inhomogeneous magnetic field acts as a source term in the evolution equations 
of the gravitational perturbations, the CMB observables may be affected. To investigate 
this situation, the relevant system of equations describing the large-scale magnetic fields is 
slightly different from the one introduced in Eqs. (|3.12|) . (|3.13jl and (|3.14|) (where, however. 
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the displacement current may be neglected). In fact, we are here interested in the dynamics 
taking place for typical frequencies much smaller than ojp^ e and ujp^ {. In this situation, 
from the two- fluid equations of Eqs. (|3.12|) . and (|3.14l) it is possible to derive various 

one-fluid plasma descriptions by defining a set of appropriate one-fluid variables such like 
the total current [21 EZ] 

J = e{niVi - UeVe), (4.1) 

the center-of-mass velocity 

V = ■ , (4.2 

me + TTii 

and so on. Among the one-fluid descriptions, a particularly simple possibility is provided by 
magnetohydrodynamics (MHD) that can be characterized by the following set of conditions: 

V-^ = 0, V-i? = 0, V-J = 0. (4.3) 

As mentioned above, the derivative of the electric field is neglected in MHD analog of Eq. 
(|3.14() . Moreover, since the current is solenoidal, as the magnetic field intensity, it follows, 
from Maxwell's equations, that 

J = X B. (4.4) 
47r 

The induced (Ohmic) electric field is therefore given, in the absence of Hall and thermo- 
electric terms, by 

E + vxB = -^VxB, (4.5) 
47r(T 

where a is the conductivity and where we used units where the speed of light is equal to 1. In 
a relativistic plasma a ~ T/ocm while in a non-relativistic plasma a ~ (T/aem)(r/^n'e)^^^- 
Since the conductivity increases as we go back in time, in the early Universe the condition 
i?-|-iTxi?/c~Ois approximately verified. Under these conditions, two (approximate) 
conservations laws may be derived, namely the magnetic flux conservation 

— f B-dT. = — f V X V X i? • dS, (4.6) 

and the magnetic helicity conservation 

d 



, , , d-^xA ■B]= -- — / d-^x B-V xB. (4.7) 
dr \Jv J 47r(T Jv 

In Eq. ()4.6() S is an arbitrary closed surface that moves with the plasma. According to 
Eq. (|4.6|) . in MHD the magnetic field has to be always solenoidal (i.e. V ■ B = 0). Thus, 
the magnetic flux conservation implies that, in the superconducting limit (i.e. a oo) 
the magnetic flux lines, closed because of the transverse nature of the field, evolve always 
glued together with the plasma element. In this approximation, as far as the magnetic field 
evolution is concerned, the plasma is a collection of (closed) flux tubes. The theorem of 
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flux conservation states then that the energetical properties of large-scale magnetic fields 
are conserved throughout the plasma evolution. 

While the flux conservation concerns the energetic properties of the magnetic flux lines, 
the magnetic helicity, i.e. Eq. ()4.7I) . concerns chiefly the topological properties of the 
magnetic flux lines. In the simplest situation, the magnetic flux lines will be closed loops 
evolving independently in the plasma and the helicity will vanish. There could be, however, 
more complicated topological situations |lU2j where a single magnetic loop is twisted (like 
some kind of Mobius stripe) or the case where the magnetic loops are connected like the 
rings of a chain: now the non-vanishing magnetic helicity measures, essentially, the number 
of links and twists in the magnetic flux lines j24j . Furthermore, in the superconducting 
limit, the helicity will not change throughout the time evolution. The conservation of the 
magnetic flux and of the magnetic helicity is a consequence of the fact that, in ideal MHD, 
the Ohmic electric field is always orthogonal both to the bulk velocity field and to the 
magnetic field. In the resistive MHD approximation this conclusion may not apply 

The quantity at the right-hand-side of Eq. ()4.7() . i.e. B ■ V x B is called magnetic 
gyrotropy and it is a gauge-invariant measure of the number of contact points in the mag- 
netic flux lines. The magnetic helicity (or the magnetic gyrotropy) and the Lorentz force 
can be used to give a general classification of the possible inhomogeneous magnetic field 
configurations that are relevant in a MHD context. Since the current density is solenoidal, 
we have that, in MHD, 

Fb = Jx B ^{V X B) X B. (4.8) 

Looking at the expressions of the magnetic gyrotropy and of the Lorentz force, there are, 
in principle, two extreme situations that may arise: 

• in the first case V x B x B = and the magnetic gyrotropy, i.e. B-V x B, is maximal; 

• in the second case the Lorentz force is maximal and the magnetic gyrotropy vanishes. 

In the case when the Lorentz force term vanishes, the configuration is said to be force- 
free. In this limit the MHD equations (coupled with the Navier-Stokes equation) become 
particularly simple and appropriate for the deduction of useful scaling laws [211 EHl IHOl 
I61j . Fully inhomogeneous magnetic field configurations may then gravitate and affect the 
evolution equations of the tensor, vector and scalar modes of the geometry. Therefore, in the 
following section, the evolution equations of metric fluctuations will be investigated in the 

^"The presence (or absence) of magnetic helicity plays an important role for the possibility of inverse 
cascades in the pre-recombination plasma. In particular, one of the evolutions of the correlation scale 
presented in Fig. Q refers to the case where the magnetic field is helical |58|. 
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presence of fully inhomogeneous magnetic fields. MHD equations in curved space-times are 
also relevant for a number of different problems arising in the physics of large-scale magnetic 
field. Concerning the various aspects of this subject see, for instance, Refs. jnHESlini] and 
references therein. 

As a last comment, we would like to stress that the magnetic helicity as well as the mag- 
netic gyrotropy are parity-odd quantities. Therefore, recalling the considerations related to 
Eq. ()3.71|) . if sources of gyrotropy and helicity are present after equality the BE and TB 
correlators may not be zero. This point has been discussed in [HSI (see, however, also |47j). 

4.1 Fully inhomogeneous magnetic correlators 

Nearly all mechanisms able to generate large scale magnetic fields imply the existence of a 

stochastic background of magnetic disturbances that could be written, in Fourier space, as 
11 

{Biik,T)B*ip,T)) = i27rfP,,{k)5^'\k-p), (4.9) 

where 

P^Jik) = Qik)(^6,, - Qik) = Qok^. (4.10) 

From Eq. 1)4. 1U() the magnetic field configuration of Eq. 1)4. 9() depends on the amplitude of 
the field Qq and on the spectral index m. 

It is easy to check that, in the case of the configuration H4.9() the magnetic gyrotropy is 
vanishing, i.e. 

{B-VxB)=0, {{B-VxBf)=0, (4.11) 

where the second expression denotes, for short, the two-point function of the magnetic 
gyrotropy. On the contrary it should also be clear that, from Eq. ()4.9() . that ((V • [V x 
B)xB])^0. 

There are situations where magnetic fields are produced in a state with non-vanishing 
gyrotropy (or helicity) [nSl EH ESI EHl • In this case, the two point function can be written 
in the same form given in Eq. (|4.9|) 

{Biik,T)B*ip,T)) = {2TifPi,{k)5^^\k-p), (4.12) 

but where now 

hj{k) = Q{k)[5,j - ^) +iQ(k)e,,ej, Qik) = Qok^. (4.13) 

^^For the Fourier transforms we use the following conventions: Bi{x) — {2tt)~"^ J d^ke^^^'^ Bi{k) and, 
conversely, B^(k,) = / cPxe^'^''^ B^{x). 



28 



From Eq. H4.13() we can appreciate that, on top of the parity-invariant contribution (already 
defined in Eqs. H4.9|l and ()4.1U|) '). there is a second term proportional to the Levi-Civita eiji. 
In Fom'ier space, the introduction of gyrotropic configurations implies also the presence of 
a second function of the momentum Q{k). 

As it will become more plausible from the considerations reported in Section 5, the 
two-point function determines, eventually, the higher-order correlation functions. Since 
the contribution of fully inhomogeneous magnetic fields to the energy-momentum tensor is 
quadratic, it is rather frequent, in various calculations, to evaluate four-point correlation 
functions. Fourier space is not always the best framework for the evaluation of higher-order 
correlation functions. Indeed, a swifter derivation of various results can be achieved, in some 
cases, by staying in real space (and by eventually Fourier transforming the final result). 
Let us therefore elaborate on the real space correlators associated with the configurations 
discussed in Eqs. and ^J^. 

The real space two-point function can be written, for instance, as 

Gijir) = {Bi{x + r)Bj{x)) = ri{r)5ij + J'2{r)rir.j + J'^{r)eijkr^ . (4.14) 

If ^3(r) = then we are in the situation described by Eq. 1)4. 9|) . If J'z{r) ^ the framework 
is the one defined by Eq. 1)4. 12(1 . 

Consider, for instance, the case T^ir) = 0. Then, using the parametrization (|4.14|) . Eqs. 
dlHj) and (jlHH) lead to 

QiAr) = 0±m, + [a|,(r) - G^irf-^, (4.15) 



where 



joikr) 



kr 



^ll(^)=(2^/^''^^^(^)^' 

where jo{kr) and ji{kr) are the usual spherical Bessel functions |69l l7L)j : 

. , sin/cr , sin kr cos kr „, 

^0(^0 = ^, nikr) = jj^-^. (4.17) 

Another useful parametrization of the real space correlators is the one where the func- 
tions Fi{r) and F2{r) are written as derivatives of a function /(r^) whose asymptotic be- 
haviour can be determined from the trace of the two-point function i 66 . More precisely, in 
the case F^(r) = 0, the functions appearing in Eq. 1)4. 14() can be writtten as 

•^iW = 7^[^'/(r^)], W = -7^[/(^')]- (4.18) 
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With the definitions ()4.18|) . Eq. H4.14() is clearly transverse (i. e. d^Gij = 0). Moreover, 
defining as Q{r) = Tr[^jj], we also have 

^^W = 2^[^'/(r')]. (4.19) 

It is often useful, in practical estimates, to regularize the two-point function by using an 
appropriate "windowing" [J^. Two popular windows are, respectively, the Gaussian and 
the top-hat functions , i.e. 

Wg(A:,L) = e-^, Wth(A:,L) = ^h{kL). (4.20) 

For instance, the regularized trace of Qij{r) with Gaussian filter can be obtained from 
the previous expressions by shifting Q{k) — > Q{k)Wg{k, L). The result is 

^" (27r)2L3+-^^ 2 '2' 4L2j' ^^'^^^ 

where -F(a, x) =i -Fi(a, 6, x) is the confluent hypergeometric function [691 170j . Notice that 
the integral appearing in the trace converges for m > —3. The amplitude of the magnetic 
power spectrum Qq can be traded for ^(0) = Bj^ (as suggested, for instance, in [72j). From 
Eq. 1)4. 21() we have that Qo becomes 

Qo = ^^^4 -Bl (4.22) 

where = 2tt/L. 

In the real space approach the higher order correlators are reduced to the calculation of 
derivatives of special functions (depending upon the regularization scheme). For instance, 
we can compute easily, in general terms the correlation function for the energy density 
fluctuations. Defining 

6pb{x) = ^B,{x)B\x), (4.23) 
it is sometimes important to evaluate 

S{r) = {5pB{x)6pB{y)) - {6pB{x)){SpB{y)). (4.24) 

By using the decomposition (|4.14l) we will have (taking for simplicity the case -F3(r) = 0), 

= ^[3F^{r) + 2r^Fi{r)F2{r)+r^Fiir) 



327r^ 
1 

32^ 



3/^(.^)+4.V(r^)^ + 2.^,^^, 
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where the second equahty follows from Eq. H4.18() . 

With these (or similar techiques) higher order correlation functions for magnetic inho- 
mogeneities may be computed for the interesting observables like the Lorentz force |72j , the 
magnetic gyrotropy [HHI and so on. 

It is appropriate to comment here on the parametrizations usually employed in the 
analysis of the various mechanisms for the magnetic field production in the early Universe. 
This topic (together with other related subjects) has been reviewed in Ref. [HS]) and, 
consequently that discussion will not be repeated here. For sake of simplicity consider the 
mechanisms relying on super-adiabatic amplification of vacuum fluctuations of the gauge 
fields (see also Fig. I^J. In this context what one is led to compute is the amplification of 
the magnetic field from some initial (vacuum) state. The precise amount of amplification 
(as well as the spectral properties of the final configuration) is determined by the features 
of the pump field that changes from model to model. Therefore, initially (i.e. for r —oo) 
the quantum mechanical operator related to the magnetic field intensity can be decomposed 
as 

^) = ^ E < / d'k KAaj;Jk{r)e-''-^^ + alj*,{T)e^'^% (4.26) 
where, within the present conventions, [a^^,a^ ] r — > +00 

Bn^. ^) = ^,T.4j d'k K[h^^^^gu{T)e~^^-' + h\^^gl{T)e^^-% (4.27) 

Since both sets of modes are complete, the old modes can be expressed in terms of the new 
ones, i.e. 

A(r) = c+(A;)5fe(r) + c.{k)gl{T). (4.28) 

Notice that /a,.(t) are normalized, as r — > —00 as e~*'^'^/ \/2A;, since we demand that the initial 
state is the vacuum. Therefore, with this normalization, C-{k), for r +00 parametrizes 
the relevant mode-mixing. 

Inserting Eq. (|4.28j) back into Eq. (|4.26|) and imposing the continuity of the operators 
we also have 

|c+(A;)|2 - |c_(A;)|2 = 1. (4.30) 
The final value of the two-point function can then be obtained as 

{omx+mm ^j^J d'k(^s,, - ^^Q{k)e-^^^-\ (4.31) 
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where now Q{k) ^ k\c^{k)\'^ . In deriving Eq. (|4.31() we used the identity 

|gfc(r)|2 = (1 + 2|c_(A:)|2)|/fc(r)|2 - c.{k)c+ik)fUrf - c*_{k)c\{k) h{Tf , (4.32) 

together with Eq. ()4.3U|) . Notice that the the second and third terms in Eq. ()4.32|) are 
typically oscillating and they have been dropped in Eq. (|4.31|) . 

Super-adiabatic amplification alone leads always to correlation functions falling in the 
general class introduced in Eq. (|4.14|) but with J--i{r) = . The primordial spectrum must 
be processed and the finite values of the thermal diffusivity as well as of the conductivity 
lead to an effective ultraviolet cut-off j63j . 

4.2 Primordial or astrophysical seeds? 

We are now are going to give some examples of magnetic fields that may be generated 
in the early stages of the life of the Universe. This topic has been already discussed in 
different review articles (see, for instance, for a short account 0Sj and [HSl for a more 
extended review). As elaborated in the introduction, the mechanisms for the generation of 
magnetic fields that may be potentially relevant for CMB physics can be divided, broadly 
speaking into two categories: mechanisms operating inside the Hubble radius (like phase 
transitions occurring, for instance, at the electroweak or at the QCD epoch) and mechanisms 
where the correlation scale of the field grows larger than the Hubble radius (as in the 
case of infiationary mechanisms). Furthermore, there is also the appealing possibility of 
so-called mixed mechanisms where the magnetic field is created, initially, in the context 
of some (conventional or unconventional) infiationary model. Then the various modes of 
the primordial spectrum of the hypercharge field will reenter at different times during the 
radiation epoch and, at this stage, will be further amplified by the processes taking place 
inside the Hubble radius. 

In the case of the electroweak phase transition (see, for instance, |62j for a comprehensive 
review on this subject) the typical strength of the produced seeds may change depending on 
the order of the phase transition (PT). We do know that, in the framework of the minimal 
standard model of particle interaction the phase transition cannot be strongly first-order 
for values of the Higgs mass larger than the mass of the W boson. This conclusion was 
reached by lattice studies of the phase diagram of the electroweak phase transition [73] . Let 
us now consider, separately, the cases when the PT is strongly first-order and second order. 

In a first-order phase transition the phases of the complex order parameter of the nu- 
cleated bubbles are not correlated. The bubbles of the broken phase expand at a velocity 
V ~ 10~^ c where c is the speed of light. At the epoch of the electroweak phase transi- 
tion the temperature of the plasma is of the order of 100 GeV while the Hubble radius. 
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i.e. ~ 3 cm. The average distance between the centers of the nucleated bubbles is 

of the order of W~^H^ |621 175j . The finite conductivity effects will give rise to diffusion 
whose consequence will be to smooth out the oscillations of the created magnetic field. The 
magnetic flux will the escape the intersection regions and penetrate the colliding bubbles 
where the evolution is governed by ordinary MHD. The magnetic fields will typically have a 
correlation scale much smaller than the Hubble radius. Even granting for a rather efficient 
diffusion of the magnetic flux, the correlation scale of the magnetic field will be at, at most, 
as large as the Hubble radius. In the absence of any form of inverse cascade we will have 
that the typical correlation scale at the onset of galactic rotation will be of the order of 10 
to 100 A.U. (recall lA.U. = 1.49 x lO^^cm). The conservative estimate of Ref. |76| is even 
smaller, i.e. the authors get a magnetic field \B\ ~ 10~^'^ G over a typical correlation scale 
of 0.1 A. U. . As suggested in the introduction, it is plausible that some form of inverse 
cascade may eventually take place between the electroweak epoch and the time of electron 
positron annihilation [SSj (see also [SSI and references therein). If the PT is of second order 
the situation is even worse. The correlation scale is likely to be of the order of 10/Tew, i-e. 
10 ^^H^^ at Te„. 

There have been also ideas concerning a possible generation of magnetic fields at the 
time of the QCD phase transition occurring roughly at T ~ 140MeV, i.e. at the moment 
when free quarks combine to form colorless hadrons. At this time the Hubble radius is of 
the order of 4 x 10^ m. The mechanism here is always related to the idea of a Biermann 
battery with thermoelectric currents developed at the QCD time Since the strange 
quark is heavier than the up and down quarks there may be the possibility that the quarks 
develop a net positive charge which is compensated by the electric charge in the leptonic 
sector. Again, invoking the dynamics of a first-order phase transition, it is argued that the 
shocks affect leptons and quarks in a different way so that electric currents are developed 
as the bubble wall moves in the quark-gluon plasma. In [77j the magnetic field has been 
estimated to be G at the time of the QCD phase transition and with typical scale of 

the order of the meter at the same epoch. 

In [7S1IZH1 it has been pointed out that, probably, the magnetic fields generated at the 
time of QCD phase transition may be much stronger than the ones estimated in (77]. The 
authors of [781 179j argue that strong magnetic fields may be generated when the broken and 
symmetric phase of the theory coexist. The magnetic fields generated at the boundaries 
between quark and hadron phases can be, according to the authors, as large as 10^ G over 
scales of the order of the meter (and possibly even larger) at the time of the QCD phase 
transition. At the onset of gravitational collapse and over a typical scale of the Mpc the 



^^See later in this Section for a general discussion of "batteries" in the framework of the generalized Ohm 
law 
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magnetic field will be of the order of 10 G with white noise spectrum (m ~ in the 
notations of this section) [75] , 

A model of variation of gauge couplings motivated by the evolution of the internal 
dimensions |^ leads to a variety of spectra that may be different depending on the number 
of extra-dimensions and on their rate of evolution (see also |821 1831 184jl. While a phase 
of expanding internal dimensions may lead to large magnetic fields, the allowed region in 
the parameter space gets reduced in the case of contracting extra-dimensions. In |SJ only 
the gauge zero- modes have been considered. In a recent calculation [S^, the role of the 
internal momenta has been taken into account with the result that also in the contracting 
case large magnetic seeds may be produced. In the case of expanding internal dimensions 
the spectral index is m ~ — 1 — rifi where n is the number of internal dimensions and /x 
is the exponent parametrizing, on conformal time, their evolution. In the case of Kasner 
models [SJ fi = \/3n/ (n + 2). If internal dimensions contract we have, effectively, m ~ 

Large-scale magnetic fields can also be produced in the context of inflationary models. 
Here the idea is that quantum fluctuations of an Abelian gauge field (for instance the 
hypercharge field) are inside the Hubble radius during inflation. Then, these fluctuations 
exit the horizon and their correlation scale grows faster than the Hubble radius itself. The 
amount of amplification achieved in this processes depends crucially on the amount of 
breaking of conformal invar iance of the evolution equation of the Abelian gauge fields. A 
lot of ideas concerning the possible breaking of conformal invariance (see, for instance, |63j 
and jHO]) have been proposed so far and it would be impossible to summarize them all in 
this topical review. 

Instead going through the various ideas it may be more useful to give further details 
about a specific way of breaking conformal invariance like the one suggested in the context 
of string cosmological models |86| . The string effective action contains a direct coupling 
between the (scalar) dilaton field and the electromagnetic field that appears even at tree 
level in the loop expansion and to lowest order in the string tension expansion. To lowest 
order the action will then be 

S = -\j d^x^ge-^Y^pY^f" , (4.33) 

where Yap is the field strength. If the gauge coupling evolves during the inflationary dynam- 
ics, then the magnetic fields will be amplified over different length-scales. In particular, in 
the model presented in [SHI, the gauge coupling, i.e. e*^/^ is a growing function of the cosmic 
time coordinate. When the Universe becomes dominated by radiation (after inflation) the 
gauge coupling freezes to a value that is of order of 10~^. More specifically, there are two 
phases that characterize the model. In the first phase the gauge coupling increases sharply 
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(the so-called dilaton driven phase). In the second phase the gauge couplmg is just slowly 
increasing (the so-called stringy phase). 

In the context of string cosmological models there are two possible ranges in the values 
of m, i.e. the spectral index of the produced magnetic fields. If the relevant modes (i.e. 
comoving scales larger than the Mpc scale) exit the Hubble radius for the first time during 
the so-called dilaton driven phase, the resulting spectrum leads to an m that could be 
written as [SS] 

m = 1= , (4.34) 

where S = /3? represents the possible effect of internal dimensions whose radii bi shrink, 
during the pre-big bang phase, like (— t)'^' for t — > 0~. In the case S = we will have 
m ~ — \/3. 

The second possibility is that the relevant large-scale modes left the Hubble radius 
during the so-called string phase. During this phase the ratio between the (cosmic) time 
derivative of the dilaton field and the Hubble rate is roughly constant. Defining as (3 the 
mentioned ratio ip/2H, we will have that m ~ 1 — 2/3. If the value of /3 gets frozen at the 
value corresponding to the last e-folds of the pre-big bang phase we will have, for instance, 
allowed values between [3 ~ 1.8 and /3 ~ 2 corresponding to m ~ —3 (i.e. —2.9, —2.8). This 
parameter can be computed once a specific model of the high-curvature (stringy) regime 
is given but it cannot be obtained, at present, from more general arguments. Notice that, 
in this case, Bl may even be of the order of 10^^*^ corresponding to a ratio between the 
magnetic and radiation energy density of the order of 10~^ |86j . 

It is also possible to generate large scale magnetic fields within the Abelian Higgs model 
where the Higgs current densities allow to obtain a (minute) amplification of the hypermag- 
netic field. The amplitude of these fields is irrelevant for applications (i.e. the produced field 
intensity is thirty orders of magnitude smaller than the most optimistic dynamo require- 
ments). The typical spectral slope is m = 1 that corresponds to the magnetic spectrum of 
the vacuum once the canonical quantum mechanical normalization is imposed on the vector 
potentials. Slightly different spectra can be obtained in the case of the models proposed in 
[571 1551 15^ where various charged scalars are considered together (see also 1^0] ) . 

It is appropriate to comment here on the possibility that large-scale magnetization is a 
purely astrophysical phenomenon. As in the case of the primordial hypothesis, it is rather 
difficult to summarize, in a reasonably self-contained manner all the hypothesis that have 
been put forward in the last fifty years. An (incomplete) list of proposals can be found, for 
instance in |631 1911 192j . The modest aim of the remaining part of this section will therefore 
be to introduce the general framework. 
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Many (if not all) the astrophysical mechanisms proposed so far are related to what is 
called, in the jargon, a battery. The idea is, in short the following. MHD, as previously 
stressed in this section, can be "derived" from a two-fluid theory by appropriately defined 
one-fluid variables that obey a set of effective Maxwell equations supplemented by the 
Ohm law (i.e. Eq. ()4.5() 1 and coupled to the Navier-Stokes equation. Now, we already 
remarked right before Eq. ()4.5() that the proposed form of the Ohmic electric field would 
contain, in principle, also a thermoelectric term that is a remnant of the subtraction of the 
momentum conservation equations that hold separately (for electrons and ions) in the two- 
fluid description from which MHD is derived. The explicit form of the generalized Ohmic 
electric field in the presence of thermoelectric corrections can be written as 

E=-,.B+^-^. (4.35) 

By comparing Eq. (|4.5|) with Eq. 1)4. 35(1 , it is clear that the additional term at the right hand 
side, receives contribution from a temperature gradient. In fact, restoring for a moment 
the Boltzmann constant ks we have that since Pe = ksUeTe, the additional term depends 
upon the gradients of the temperature, hence the name thermoelectric. It is interesting 
to see under which conditions the curl of the electric field receives contribution from the 
thermoelectric effect. Taking the curl of both sides of Eq. 1)4. 35() we obtain 

VxE=- — V^B + V{vxB) - = -^, 4.36 

where the second equality is a consequence of Maxwell's equations. From Eq. (|4.36|) it is 
clear that the evolution of the magnetic field inherits a source term iff (i.e. if and only if) 
the gradients in the pressure and electron density are not parallel. If VPc 1 1 Vric a fully 
valid solution of Eq. (|4.36|) is i? = 0. In the opposite case a seed magnetic field is naturally 
provided by the thermoelectric term. The usual (and rather general) observation that one 
can make in connection with the geometrical properties of the thermoelectric term is that 
cosmic ionization fronts may play an important role. For instance, when quasars emit 
ultraviolet photons, cosmic ionization fronts are produced. Then the intergalactic medium 
may be ionized. It should also be recalled, however, that the temperature gradients are 
usually normal to the ionization front. In spite of this, it is also plausible to think that 
density gradients can arise in arbitrary directions due to the stochastic nature of density 
fluctuations. 



^^For simplicity, we shall neglect the Hall contribution arising in the generalized Ohm law. The Hall 
contribution would produce, in Eq. I|4.35|l a term J x B/rice that is of higher order in the magnetic field and 
that is proportional to the Lorentz force. The Hall term will play no role in the subsequent considerations. 
However, it should be borne in mind that the Hall contribution may be rather interesting in connection with 
the presence of strong magnetic fields like the ones of neutron stars (i.e. 10^"^ G). This occurrence is even 
more interesting since in the outer regions of neutron stars strong density gradients are expected. 
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In one way or in another, astrophysical mechanisms for the generation of magnetic fields 
use an incarnation of the thermoelectric effect. In the sixties and seventies, for instance, it 
was rather popular to think that the correct "geometrical" properties of the thermoelectric 
term may be provided by a large-scale vorticity. As it will also be discussed later, this 
assumption seems to be, at least naively, in contradiction with the formulation of inflationary 
models whose prediction would actually be that the large-scale vector modes are completely 
washed-out by the expansion of the Universe. Indeed, all along the eighties and nineties 
the idea of primordial vorticity received just a minor attention. 

The attention then focused on the possibility that objects of rather small size may 
provide intense seeds. After all we do know that these objects may exist. For instance the 
Crab nebula has a typical size of a roughly 1 pc and a magnetic field that is a fraction of the 
m G. These seeds will then combine and diffuse leading, ultimately, to a weaker seed but 
with large correlation scale. This aspect, may be, physically, a bit controversial since we 
do observe magnetic fields in galaxies and clusters that are ordered over very large length 
scales. It would then seem necessary that the seed fields produced in a small object (or in 
several small objects) undergo some type of dynamical self-organization whose final effect 
is a seed coherent over length-scales 4 or 5 orders of magnitude larger than the correlation 
scale of the original battery. 

An interesting idea could be that qualitatively different batteries lead to some type of 
conspiracy that may produce a strong large scale seed. In PJ21 it has been suggested that 
Population III stars may become magnetized thanks to a battery operating at stellar scale. 
Then if these stars would explode as supernovae (or if they would eject a magnetized stellar 
wind) the pre-galactic environment may be magnetized and the remnants of the process 
incorporated in the galactic disc. In a complementary perspective, a similar chain of events 
may take place over a different physical scale. A battery could arise, in fact in active galactic 
nuclei at high red-shift. Then the magnetic field could be ejected leading to intense fields in 
the lobes of "young" radio-galaxies. These fields will be somehow inherited by the "older" 
disc galaxies and the final seed field may be, according to j92j as large as 10~^ G at the 
pre-galactic stage. 

Up to now only a necessarily reduced account of the primordial and of the astrophysical 
hypotheses has been presented. In spite of the mechanism leading to a seed field, the crucial 
problem will be how to go from the seed to the actual galactic (or even cluster) magnetic 
field. The galactic rotation period is of the order of 3 x 10^ yrs. This scale should be 
compared with the typical age of the galaxy. All along this rather large dynamical time-scale 
the effort has been directed, from the fifties, to the justification that a substantial portion 
of the kinetic energy of the system (provided by the differential rotation) may be converted 



37 



into magnetic energy amplifying, in this way, the seed field up to the observed value of the 
magnetic field, for instance in galaxies and in clusters. In recent years a lot of progress has 
been made both in the context of the small and large-scale dynamos. This progress was 
also driven by the higher resolution of the numerical simulations and by the improvement in 
the understanding of the largest magnetized system that is rather close to us, i.e. the sun. 
More complete accounts of this progress can be found in the second paper of Ref. and, 
more comprehensively, in Ref. j94j . The main ingredients of the large-scale dynamo can be 
also found in Apart from the aspects involving solar physics and numerical analysis, 
better physical understanding of the role of the magnetic helicity in the dynamo action has 
been reached. This point is crucially connected with the two conservation laws arising in 
MHD, i.e. the magnetic flux and magnetic helicity conservations whose relevance has been 
already emphasized, respectively, in Eqs. 1)4. 6|1 and 1)4. 7() . Even if the rich interplay between 
small and large scale dynamos is rather important, let us focus on the problem of large-scale 
dynamo action that is, at least superficially, more central for the considerations developed 
in the present paper. In short, the picture of large-scale dynamos can be summarized as 
follows. In a turbulent environment the MHD equations can be appropriately averaged 
with the result that the mean magnetic seed is amplified exponentially fast. This approach 
is described in detail in [5^ (see also USD- The most important point to bear in mind 
is that the turbulent velocity field (whose correlation scale is typically shorter than the 
correlation scale of the magnetic seed) must lack mirror symmetry. In other words the 
velocity field must be parity-odd. This occurrence necessarily implies that the (averaged) 
magnetic diffusivity equation inherits a term that is proportional to the curl of the magnetic 
field though a coefficient, the celebrated a term, that would vanish if the velocity field would 
be parity-even. It can be demonstrated that the a term is indeed proportional, in turn, to 
the average kinetic helicity, i.e. {v ■ V x v). The a term decreases more slowly, at large 
length-scales, than the diffusivity terms. Therefore, for sufficiently large scales the magnetic 
seed can be efficiently (exponentially?) amplified. For long time this oversimplified picture 
was thought to be essentially consistent even if computationally challenging. 

Already at a qualitative level it is clear that there is a clash between the absence of 
mirror-symmetry of the plasma, the quasi-exponential amplification of the seed and the 
conservation of magnetic flux and helicity in the high (or more precisely infinite) conduc- 
tivity limit. The easiest clash to understand, intuitively, is the flux conservation versus the 
exponential ampliflcation: both flux freezing and exponential ampliflcation have to take 
place in the same superconductive (i.e. 0) limit. The clash between helicity con- 

servation and dynamo action can be also understood in general terms: the dynamo action 
implies a topology change of the conflguration since the magnetic flux lines cross each other 
constantly [55] . 
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One of the recent progress in this framework is a more consistent formulation of the 
large-scale dynamo Droblem |931l94j : large scale dynamos produces small scale helical fields 
that quench (i.e. prematurely saturate) the a effect. In other words, the conservation of 
the magnetic helicity can be seen, according to the recent view, as a fundamental constraint 
on the dynamo action. In connection with the last point, it should be mentioned that, in 
the past, a rather different argument was suggested it was argued that the dynamo 
action not only leads to the amplification of the large-scale field but also of the random field 
component. The random field would then suppress strongly the dynamo action. According 
to the considerations based on the conservation of the magnetic helicity this argument seems 
to be incorrect since the increase of the random component would also entail and increase 
of the rate of the topology change, i.e. a magnetic helicity non-conservation. 

In summary we can therefore say that: 

• both the primordial and the astrophysical hypothesis for the origin of the seeds de- 
mand an efficient (large-scale) dynamo action; 

• due to the constraints arising from the conservation of magnetic helicity and magnetic 
flux the values of the required seed fields may turn out to be larger than previously 
thought at least in the case when the amplification is only driven by a large-scale 
dynamo action 

• magnetic flux conservation during gravitational collapse of the protogalaxy may in- 
crease, by compressional amplification, the initial seed of even 4 orders of magnitude; 

• compressional amplification, as well as large-scale dynamo, are much less effective in 
clusters: therefore, the magnetic field of clusters is probably connected to the specific 
way the dynamo saturates, and, in this sense, harder to predict from a specific value 
of the initial seed. 



5 Magnetized metric perturbations 

Consider a conformally flat background geometry 

9p.v = (5.1) 

where rj^i, = diag(l, —1, —1, —1) is the flat Minkowski metric. Then the fluctuations of the 
metric may be written as 

^Ofiu = Stgfiu + SvOnu + SsOfiu, (5.2) 
^*The situation may change if the magnetic fields originate from the combined action of small and large 
scale dynamos like in the case of the two-step process described in 1921 . 
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where 5t, 5^ and 5s denote the tensor, vector and scalar fluctuation of the geometry. 

The evolution of the various modes of the geometry are sourced by the corresponding 
modes of the energy-momentum tensor of the fluid sources and of the magnetic field. In 
particular, denoting with T^^, the energy-momentum tensor of the magnetic field we will 
have, following the notation of Eq. 1)5. 2|1 . that 

ST^iy = 6tTf^u + Sv'T',_iu + SsTf^u- (5.3) 

The fluctuations ST^i, are quadratic in the intensity of the fully inhomogeneous magnetic 
flelds. This means, that a non-Gaussian signal should be, in general, expected. The rationale 
for this statement is that while the fluctuations of fully inhomogeneous magnetic fields may 
be Gaussian (see, for instance, Eq. 1)4. yi) ). the fluctuations of the energy- momentum tensor 
will be distributed as the square of a Gaussian variable and will therefore be non-Gaussian. 
Recently, more accurate calculations of the possible non-Gaussian effects arising in the 
presence of fully inhomogeneous magnetic flelds have been presented |97| . 

In the following the evolution equations of tensor, vector and scalar modes will be 
treated within a (first-order) "Bardeen" approach [52] (see also |4U| for an introduction to 
the key features of the problem in the absence of magnetic fields). The evolution of fully 
inhomogeneous magnetic fields can be usefully treated also within a covariant formalism. 
For this type of analysis see the interesting studies reported in Refs. [Ml IIUUI llUlj . 



5.1 Magnetized tensor modes 

The tensor fluctuations of the metric deflned as 

5tgij = -a^h,j, dih] = 0, l4 = 0, (5.4) 

are automatically invariant under infinitesimal coordinate transformations. Equation 1)5. 4() 
implies that hij carries two independent degrees of freedom. From the fiuctuation of the 
Ricci tensor it is easy to obtain that the evolution equations are 

h{" + 2nh( - V^h{ = -16TTGa^5tT->, (5.5) 

where, as mentioned before, 5tT-'' is the transverse and traceless component of the energy- 
momentum tensor of the electromagnetic fields. Given the specific magnetic field config- 
uration, Eq. ()5.5|) allows to determine the amplitude of the gravitational waves. The 
determination of hij allows the determination of the tensor Sachs- Wolfe contribution, i.e. 

h'ijn'n^dr. (5.6) 









2/, 
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In the tensor case the fluctuations of the geometry do not mix with the fluctuations of the 
fluid sources. Furthermore, in the tensor case, as it is well known, only the Sachs- Wolfe 
integral contributes. Tensor modes have been discussed in a number of papers. In and 
jlL)2j the possible stochastic backgrounds of gravitational radiation stemming from inho- 
mogeneous magnetic fields have been computed with particular attention to configurations 
carrying non- vanishing magnetic gyrotropy. Recently gravitational waves produced through 
a phase of helical turbulence have been discussed in |lU3j . For a discussion of the interplay 
of gravitational waves and inhomogeneous magnetic fields see also |1U4| . 

5.2 Magnetized vector modes 

The evolution equations arising in the case of vector modes are a bit more complicated than 
in the case of tensor modes. The vector modes of the geometry are parametrized in the 
following way 

5v5o* = -a^ft, <5v5*i = a\diWj + djWi), (5.7) 
where SjQ* = and diW^ = 0. 

The perturbed Einstein and covariant conservation equations read 

d^R"^ = 87rG(5vr; + d^rp, (5.8) 

V^^vT^" + V^^vT^'^ = 0, (5.9) 

where 

<5vr; = ip + p){u^5,v!' + S^u^u"), (5.10) 

-^vT; = ^ (-F^„F- + . (5.11) 

Equations (|5.8|) and (|5.11jl have to be supplemented by Maxwell's equations together with 
the generalized Ohm's law 

f = a^Ff^'u, (5.12) 



Recall first that 



where 



2a^ 

= :^[{diV^ + d^Vi)' + 2nidiV^ + d^V,)]. (5.13) 



V = Q + W'. (5.14) 
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The vector V is gauge-invariant. In fact, for infinitesimal gauge transformations preserving 
the vector nature of the fluctuation we have, indeed. 



Qi ^ Qi — Qi Cii 

Wi^Wi = Wi + Ci. 



(5.15) 
(5.16) 
(5.17) 



From Eq. (|5.10() . recalling that 



uo6yu' = -iV + Q'), 



we simply have 



6,T^ = -ip + p){V + Q^). 
Equation ()5.11|) become, in explicit terms. 



where the following conventions have been adopted 



Fqi — a £i, F' 



Fij — o €ij]fB , F 
E = c?£, B = a^B. 



7Z 



(5.18) 
(5.19) 

(5.20) 

(5.21) 



(5.22) 



Therefore, the perturbed set of equations can be written, in the gauge Wi = 0, 



V\Q> + 2rLQ) = ^^Fb{x) 



V 



m + 



where 



p' + p' 

P + P . 



1 



V + 



Fb{x) 



0, 



Fb = —(V X B) xB. 



(5.23) 
(5.24) 

(5.25) 
(5.26) 



Notice that the divergence- full part of the Lorentz force (i.e. V • Fb) does not contribute 
to the evolution equations of the vectors but it does contribute to the evolution equations 
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of the scalar modes. The second term at the right hand side of Eq. ()5.23|) contains, in the 
denominator, a = aac, i.e. the rescaled conductivity. This term, coming from the MHD 
form of the Poynting vector, is neghgible in the infinite conductivity hmit and it is usually 
dropped. 

Formally, the solution of the above system of equations can be written as 

Finally, the vector contribution to the Sachs- Wolfe effect can be written, in the gauge Wi = 
as 

/AT\ ^ 1 /"^/ 

[-jrj =[-V-n];( + -j {diQj+d,Qi)n'n^dT. (5.28) 

In the absence of large-scale magnetic field, the only solution of the above system is the 
one denoted by the integration constant C. In this case, if the Universe expands, the vector 
modes can only decay (see, however, |1U51 11U6T llU7p . 

In a series of papers Barrow and Subramanian jlU81 llU9j analyzed the effects of tangled 
magnetic fields both on the temperature and polarization anisotropics. More recently the 
vector modes have been analyzed numerically in Ref. jllUj (see also The obtained 

results are in qualitative agreement with the semi-analytical estimates presented in jlU81 
llU9j . In |llUj the effect of the neutrinos has been accurately parametrized in the evolution 
equations. In the equations for the vector modes written above, the (divergenceless) velocity 
field V has been taken generic. However, when the plasma content before decoupling is 
specified, V, i.e. the total velocity field will receive contribution from all the species present 
in the plasma. A similar situation will be encountered in the case of scalar modes and will 
be briefly discussed there. 



5.3 Magnetized scalar modes 

The evolution of scalar fluctuations is, technically, more complicated than the vector and 
tensor problems. In the case of scalar fluctuations, the magnetic fields enter the evolution 
equations through the magnetic energy density and through the divergence of the Lorentz 
force. Furthermore, in the scalar case, the possible gauge choices increase with respect to 
the case of the vector modes. 

In the absence of magnetic fields, we do know that the most general initial conditions for 
CMB anisotropics consist of one adiabatic mode and four non-adiabatic (or isocurvature) 
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modes (see, for instance, HOI)- In the case of adiabatic modes the fluctuations in the specific 
entropy vanish for typical length-scales larger than the Hubble radius. This occurrence 
implies a specific relation among the density contrasts of the various species present in the 
plasma. Consider, indeed, the simple case of a system formed by cold dark matter (CDM) 
and radiation. In this case, the specific entropy ? and its relative fluctuations S can be 
written as 

, S = - = ^5^-5c. (5.29) 

where ric is the number density of CDM particles; S-y and 5c are the density contrasts of the 
fluctuations in radiation and in CDM. If 5 ~ for scales larger than the Hubble radius, 
the fluctuations are said to be adiabatic. 

This naive definition can be generalized in various ways. One can, for instance, give a 
more general gauge-invaraint definition by introducing two generic fluids, the a-fluid and 
the b-fluid. In this case it is rather natural to deflne the two gauge-invariant variables 

Ca = -V'-H%, (5.30) 

Pa 

Cb = -V'-W^, (5.31) 

Ph 

that can be interpreted, respectively, as the curvature fluctuations on uniform a-density 
and b-density hypersurfaces (see, for instance, |112j and references therein). The (total) 
curvature fluctuations on uniform density hypersurfaces are then 

C = 4Ca + 4Cb^-V'-W^, (5.32) 
P P P 

where /o = Pa + Pb. From Eqs. 1)5. 3U() and ()5.31|) it is possible to define the generalized 
entropy fiuctuations as 

S = -3(Cb - Ca) = rr^ - -r^, (5.33) 

where the second equality follows from the first equality by using the covariant conservation 
equation for the each of the two fluids; w^, and are, respectively, the barotropic indices 
for the a-fluid and the b-fluid. In the case Wh = w-y = 1/3 and Wa = Eq. ()5.33|) reproduces 

Eq. dsisni). 

If we now consider the plasma after equality but before decoupling, the total entropy 
fluctuations will vanish iff 

5j^5,^ ^5c ~ ^<5b (5.34) 

where the density contrast in the neutrinos and in the baryons, i.e. 5^ and 5^, have been 
introduced. 
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If magnetic fields are included in the game, then the adiabatic mode, together with the 
four isocurvature modes, can be generalized. The logic will then be the following: 



• the evolution equations of the whole system must be written including the contribution 
of the large-scale magnetic fields; 

• the system includes, on top of the magnetic field and of the scalar fluctuations of the 
geometry also the neutrinos, the photons, the baryons and the CDM particles; 

• the equations must be solved when the relevant modes are larger than the Hubble 



• the obtained solutions will serve as initial conditions for the lowest multipoles of the 
Boltzmann hierarchy. 

It is difficult, in the scalar case, to go trough all the previous steps without changing 
gauge. Indeed it is well known, also in the absence of large-scale magnetic fields, that some 
isocurvature modes (for which the presence of anisotropic stresses is crucial) are regular in 
the synchronous gauge but are divergent (at early times) in the longitudinal gauge. This 
is just an example and, in the following, the magnetized adiabatic mode will be mainly 
discussed, quoting, when appropriate, some results of more technical investigations ||113j . 

5.3.1 Magnetized adiabatic mode 

Consider, for simplicity, the system of scalar fiuctuations in the longitudinal gauge where 
the metric fluctuations can be written in terms of the two functions (j) and ip, i.e. 



In the following, the relevant equations of the system will be discussed starting with the 
evolution equations of the fluctuations. As mentioned before, the plasma content is assumed 
to be formed by neutrinos, photons, baryons, CDM particles and fully inhomogeneous 
maggnetic fields. 

In the gauge of Eq. (|5.35() . the Hamiltonian constraint reads 



radius; 



(5.35) 



(5.36) 



where, for Ny species of massless neutrinos 




1 + 



R 




(5.37) 
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so that Ru and i?^ represent the fractional contributions of photons and neutrinos to the 
total density at early times deep within the radiation-dominated epoch. 

In Eq. (|5.36|1 the contribution of the magnetic energy density has been parametrized, 
in Fourier space, as 

nBik,v) = ^ = / d^P B^{\p- k\)B\p). (5.38) 

The appearance of convolutions in the expressions of the magnetic energy density is a 
direct consequence of the absence of a uniform component of the magnetic field whose 
background contribution, as repeatedly stressed, is vanishing. Notice also that p appearing 
in Eq. 1)5. 38(1 is the energy density of the dominant component of the fluid, so, for instance, 
in the radiation-dominated epoch p = p^, a^p^ ~ constant and Q.b — constant. 

While the Hamiltonian constraint relates the metric fluctuations (and their first time 
derivatives) to the density contrasts of the various species, the momentum constraints relates 
the metric fiuctuations (and their first derivatives) to the divergence-full part of the peculiar 
velocity field and to the MHD expression of the Poynting vector^^: 



^{R,e, + (1 - R,)e^) + + 0bf^b + 



(5.39) 



where J--Q{k,rj) is the Fourier transform of the generalized Lorentz force, i.e. 

V • [(V X B) X S] = j d^kJ^B{ky^-^ , (5.40) 
and J-'sik) is given by the following convolution 

I d^p[{k-p)B^{p)B\\k-f^)-{h-pi)B0-f^)B,[\k-m\p)]- (5.41) 
Notice that a consistent treatment of fully inhomogeneous magnetic fields implies that 

in order not to over-close the Universe. In the ideal MHD limit, i.e. a ^ oo, the contribution 
of J-'B{k) disappears from Eq. (|5.39|) . 

Another important condition on the difference of the conformally Newtonian potentials 
comes from the (i 7^ j) components of the perturbed Einstein equations: 



k'{<i)-^) = In^ 



4 

-^(yuRu - o"B(fe) 



(5.43) 



^In the following 0^ = ftul, 9c = divl and so on for the two remaining species, i.e. 6\, and 9y. 
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having defined 

The quantity ctb can be interpreted as the anisotropic stress arising thanks to the inho- 
mogeneous magnetic field and, in the force-free case, ctb = ^^b/S. Equation H5.43() relates 
the difference of the two metric fluctuations to a combination of the neutrino and magnetic 
anisotropic stresses. 

Let us finally come to the equation stemming from the perturbed {i = j) components 
of the Einstein equations: 

1.2 -, 

V'" + {2^p' + (p')n + {2n' + n^)(i) -—{(t)-^) = -n^[{Rj^ + (i - Ru))s^ + nsik)], (5.45) 

where the contribution of the other species vanishes since both baryons and CDM particles 
have a vanishing barotropic index, i.e. wy^ = Wc = 0. 

Equations (|5.36j) . (|5.39|) . (|5.43|) and (|5.45j) should be supplemented by the evolution 
equations of the various plasma quantities. For doing so, there are two complementary 
approaches. In the first approach, one may notice that the perturbed Einstein equations 
only involve the density contrasts and the peculiar velocities of the various species, i.e. the 
monopoles and the dipoles of the (perturbed) phase spece distributions. Therefore, the 
simplest approach is to describe the evolution of the different species of the plasma through 
an appropriate set of fluid equations (following from the covariant conservation equations 
for each species). The exception in this statement is represented by the neutrinos whose 
evolution includes also the quadrupole of the neutrino phase space distribution, i.e. a,^. 
To treat neutrinos one should include, therefore, also higher multipoles in the Boltzmann 
hierarchy. Alternatively, one may treat all the species through a Boltzmann hierarchy. This 
is the second, complementary, approach to the problem of initial conditions. 

Proceeding along the improved fluid approach the covariant conservation equations for 
photons and baryons imply, respectively, 

e{, + ne^ = + + ^9laneXea^{0^ - ^b) = (5.46) 

= 3V' - ^b- (5.47) 

and 



5'^ = - e^ + A^;', (5.48) 



1.2 

g'^ = —§^ + k'^(j)-\-ane(JT{Ob-e^). (5.49) 

In Eqs. ()5.46|) and ()5.49|) Thompson scattering term has been added. At early times, 
baryons and photons are synchronized so well that, to lowest order in the tight coupling 
approximation, (Tt ^ oo and 9-., — ^b- 
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The evolution of the CDM particles follows from the fluctuations of the covariant con- 
servation equation whose first-order fluctuation leads, in Fourier space, to 



6', = 3iP' - 6,. 



(5.50) 
(5.51) 



To describe neutrinos during and after horizon crossing requires a Boltzmann hierarchy 
for 0u and for the higher multipole moments, i.e. ^ > 2, of the neutrino phase-space 
density J^^e- With this caveat in mind, after neutrino decoupling, at temperatures of about 
1 MeV, massless neutrinos obey, in Fourier space, the following set of equations 



6u 



5:, = --e, + 4i{;', (5.52) 

(5.53) 

15- 10-™- <'^-''" 
where = ^^^2/2 is the quadrupole moment of the (perturbed) neutrino phase-space 
distribution and, as introduced above, J^^e is the i-th multipole. 

The closed system formed by the perturbed Einstein equations and by the fluid equations 
may be solved. For illustrative purposes, consider the magnetized adiabatic mode in the 
force-free case. The solution can be obtained as an expansion in kr which is small, i.e. 
kr < 1 when the relevant modes are outside the horizon. In the presence of a magnetic field 
the density contrasts turn out to be 



■>h,c 



>h,c 



69 - 61 -R - 8 fig 
60 {25 + 2R^) 
237-H52fi^ + 16 Rl ' 

90 {25 + 2R^) 
375 - 207 R^ - 152 Rl - 16 fi^^ 



nuk^T\ 



^Bk^T 



2^2 



(5.55) 
(5.56) 
(5.57) 



90 Ru {2b + 2Ru) 

where the barred quantities denote the adiabatic mode in the absence of magnetic field. 
In particular, for the density contrasts in the absence of inhomogeneous magnetic field we 
have 



'Su = -2(f>o 



525 + 188 i?^ + 16ii2) 

-^(f)ok^T^, 



60 (25 + 2i?^) 
(525 + 188 i?^ + 16 Rl 



(5.58) 
(5.59) 



45 (25 + 2fi^) 

where 4>o is the constant value of (p to lowest order in kr. While in the absence of magnetic 
field the adiabaticity condition of Eq. (|5.34j) is verified both to lowest and to higher order in 
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kr, in the presence of magnetic field the adiabaticity condition is violated to next-to-leading 
order. 



With similar notations we can also report the evolution for the metric fluctuations 

'6-8R^ + 2Rl 



ip = ip 



45 {25 + 2R^) 
69 - 61 i?^ - 8 Rl 



17f 



7 2 2 
k T , 



180 {25 + 2R^ 

where, as before, (/> and ^/J denote the standard adiabatic solution, i.e. 



(75 +U R^-8 Rl) 
90 {2h + 2Ru) 



V2 



^ {7b + 79R, + 8Rl) ^ 2 2 



(5.60) 
(5.61) 

(5.62) 
(5.63) 



90 {25 + 2R^) 

As in the standard adiabatic case, ipQ is determined in terms of i^o by the neutrino anisotropic 
stress. In the magnetized case the relation is modified as 



Tpo = (/)o -h -R^ 
The solution for the velocity fields is finally 



iRu-1). 



4 40 (25 + 2i?) ^ 



OBi?7^2 

4 R^ 



45 - 29 i?^ 



72 R^ (25 + 2R^ 



kW 



Oc = 0h 



6-8Ru + 2Rl 



k'r^ 



.180 {25 + 2R^) 

where, following the same notations used before, the standard adiabatic solution is 



(5.64) 

(5.65) 
(5.66) 
(5.67) 

(5.68) 
(5.69) 
(5.70) 

20 (25 + 2i?^)""-^ ■ ' ^^""^^^ 

Notice that while to leading order all the peculiar velocities of the plasma coincide, to 
next-to-leading order they are different: this occurs already in the standard adiabatic case. 

Since the adiabaticity condition is only realized to lowest order in kr it is legitimate 
to name the solution presented in Eqs. (|5.55j) - (|5.67jl quasi- adiabatic. The corrections to 



= ^k\ 



^k\ 



2 36 {25 + 2 R^y ' 

% ,2 (75 + 14 i?, - 8 Rl) 4 3 



360 {25 + 2R^) 
(75 + 14 i?^ - 8 Rl) 
360 (25 + 2 R^) 
(25 + 8i?, 



V, 
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adiabaticity are of order r^B^^r^, so they are suppressed outside the horizon and also, by 
virtue of Eq. ()5.42|) . by JIb- Notice that three regimes emerge naturally. The quasi- 
adiabatic regime where 0,B{k) < 4>o{k), the isocurvature regime ^B{k) > (poik) and the 
fully adiabatic regime, i.e. J7b(^) ^ 0o(^)- 

It is useful to recall that the presence of fully inhomogeneous magnetic fields changes 
the qualitative features of the tight coupling expansion [113j . As it is well known, the idea 
of the tight-coupling approximation is to tailor a systematic expansion in powers of the 
inverse optical depth, i.e. 1/e' in the notations employed Eq. 1)3.40(1 . Instead of being 
interested in the brightness perturbations related to Q and U we shall be interested in the 
brightness perturbation related to I, i.e. the intensity of the scattered radiation. To derive 
the evolution equations in the tight-coupling approximation for a magnetoactive plasma we 
recall that the equation for Aj can be written as 



A[ + ikfiAi = (ip' - ikiicj)) + r' 



-Ai + Aio + Aifb - ]^P2{^J)SQl 



(5.72) 



where we defined fb = 6\^/{ik) and where 5q is the same quantity already introduced in 
Eq. (|3.51() . The equation for v\y can be readily obtained from Eq. (|5.46() and it is given by 



t' 



vi + Hvb = -ik(j) [3iAii + Vb] - iWB(k), (5.73) 

a 



where 



= 4^' ^'-''^ 
a = ^^. (5.75) 

4/07 

It should be noticed that even if the magnetic field does not appear directly in Eq. (|5.72j) . it 
does appear in Eq. H5.73() . This implies that the equation for Ai (already to zeroth order in 
the tight-coupling expansion) inherits an extra source term that depend on the divergence of 
the Lorentz force. This statement can be verified by expanding, consistently, the evolution 
equations for the brightness perturbations in powers of l/|e'| |113j . For instance, to zeroth 
order in the tight-coupling expansion the (decoupled) evolution equation for the monopole 
can be writtten as 



-// , a' k"^ 
a + l ^° ^ 3(a + 1) 



Aio + — ^Aio + , .. Aio 



a + 1 3 



a J^B 
a + l 127rpba'^ ' 



(5.76) 



In the limit J^b — > the usual decoupled equation for the evolution of the monopole is 
clearly recovered. 
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5.4 Constraints on fully inhomogeneous fields 



In the following we will refer exclusively to the fully inhomogeneous case since the uniform 
case has been already discussed in Section 3. In |115| the effects of a fully inhomogeneous 
magnetic field has been considered semi-analytically. The analysis of |115j was motivated 
by the results of where, in a specific model, quite sizable large-scale magnetic fields could 
be onbtained. In particular, the energy spectrum could be, in that case, quasi-fiat, i.e. in 
the notations of Section 4, m ~ —2.9. It was then noted that a rough limit on the magnetic 
field intensity could be obtained, i.e. < 10~^ where we follow the notations of Section 
4 and take as comoving scale of regularization the Mpc. 

A semi-analytical analysis of the tensor and vector modes performed in [72]. The authors 
discussed correctly two opposite regime. In the case m < —3/2 the two-point function of 
the Lorentz force is dominated by the large-scale features of the magnetic field. On the 
contrary, if m > —3/2 the two-point function of the Lorentz force is dominated by the high 
A:-modes and hence by diffusion. In fact, the power-spectrum of the Lorentz force term can 
be estimated as 

2 



P{k) ~ 



7r(2m + 3) 



(27r)"^+35| 



2m+3 / t,\3 ^ / h \ 2n+3n 

, (5.77) 



k^J [k^J m + 3[k[) 



where k^ is the diffusion scale. The diffusion scale can be estimated, for instance, from the 
considerations of Ref. ^HI where it was argued that the diffusion length-scale is basically 
determined by the Alfven velocity times the Silk length-scale. Prom Eq. ()5.77|) we can 
see that the second term dominates for —3 < m < —3/2 while the first term dominates for 
m > —3/2. For m = —3/2 both terms may be of the same order. Notice that, in Eq. ()5.77|) . 
k^ = 2'k/L and L is the regularization length-scale introduced at the end of the previous 
Section (see Eq. ^^). 

In |72] the TT, EE, TE and BB correlations for vector and tensor modes have been 
estimated semi-analytically. Using 10% estimates the authors of Ref. 72_ obtain for m ~ —3 
(e.g. m ~ —2.9) a limit Bl < 7 x 10^^ G. This type of limit becomes more stringent as m 
increases: for m = we have < 10"^° G, while for m ~ 2 we would have < 10"^^ G. 

Recently in jlllj the calculation of the vector and tensor modes has been pushed to 
d. ~ 2000. It is clear from the reported results jlllj that for a magnetic field with Bl ~ 
10~^ G and m ~ —2.9 the TT correlations induced by the vector modes become of the same 
order to the TT power spectrum (in the absence of magnetic field) for i > 2000. Expected 
numerical disagreements with the semi-analytical estimates of |72| have been also found for 
(. > 100. In Fig. IS] the results obtained in the case of vectors for a typical magnetic field 
^^See the discussion in Section 3 on the cross power spectra whose general definition is given in Eq. Il3.fcifcil . 
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10' 10' 10= 



Figure 5: CMB power spectra (vector modes) for Bl ~ 3 X lO^^G. Adapted from Ref. 
jllUj . TT correlations in the absence of magnetic field (curve (a)); TT correlations in the 
presence of magnetic fields (curve (b)); absolute value of TE correlation (curve (c)); BB 
correlations (curve (d)); EE correlation (indicated by the arrow). These correlations have 
been defined in Section 3, see, in particular, the discussion from Eq. H3.61() to Eq. ()3.71|) : 
the notation EE, BB, TE, and so on, refers to the cross power spectra defined in Eq. (|3.66() 
(see also the text). 
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Bl ^ 3 X IQ-^G are reported (see HTU]). 

Scalar modes have been less studied than the vector and tensor modes. The expectation 
is that scalar modes do not give stronger constraints than vectors or tensors. It would 
be, however, desirable to bring the analysis of the scalar problem to the same standard of 
ordinary CMB calculations where the adiabatic mode and the different isocurvature modes 
are constrained to a given confidence level. For instance the results reported by |114j are 
difficult to intepret since bounds on the magnetic field intensity are claimed but the initial 
conditions for the analysis are not clearly specified. 

5.5 Other effects 

On top on the effects of large-scale magnetic fields on anisotropies there are effects directly 
connected with CMB spectrum like the possible distorsions. These have been studied in 
Ref. jll7j (see also jll8j l. The possible distortions of the Planckian spectrum of CMB are 
usually discussed in terms of a chemical potential which is bounded, by experimental data, 
to me < 9 X 10~^. Magnetic field dissipation at high red-shift implies the presence of a 
chemical potential. Hence bounds on the distorsion of the Planckian spectrum of the CMB 
can be turned into bounds on the magnetic field strength at various scales. In particular 
jll7j the obtained bound are such that B < 3x 10"^ G for comoving coherence lengths 
between 0.4 kpc and 500 kpc. 

As in the case of synchrotron emission also Faraday rotation measurements can be 
used as a diagnostic for foreground contamination. The idea would be, in this context to 
look for cross-correlations in the Faraday rotation measure of extra-galactic sources and 
the measured microwave signal at the same angular position. A recent analysis has been 
recently reported [IHI (see also fT^IHTlIT^ ). 

6 Concluding remarks 

Simple logic dictates that if the origin of the large-scale magnetic fields is primordial (as 
opposed to astrophysical) it is plausible to expect the presence of magnetic fields in the 
primeval plasma also before the decoupling of radiation from matter. As it was argued in 
the present topical review, CMB anisotropies are germane to several aspect of large-scale 
magnetization. CMB physics may be the tool that will finally enable us either to confirm or 
to rule out the primordial nature of galactic and clusters magnetic fields seeds. In the next 
five to ten years the forthcoming CMB precision polarization experiments will be sensitive 
in, various frequency channels between abovelOO GHz and, roughly below 900 GHz. The 
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observations will be conducted both via satellites (like the Planck satellite) and via ground 
based detectors (like in the case of the QUIET arrays). In a complementary view, the SKA 
telescope will provide full-sky surveys of Faraday rotation that may even get close to 20 
GHz. 

In an optimistic perspective the forthcoming experimental data together with the steady 
progress in the understanding of the dynamo theory will hopefully explain the rationale 
for the ubiquitous nature of large-scale magnetization. In a pessimistic perspective, the 
primordial nature of magnetic seeds will neither be confirmed nor ruled out. The plausibility 
of (highly speculative) theoretical scenarios for the origin of primordial magnetic fields 
cannot be decided only using the fertile imagination of theoretical physicists. It would 
be wiser, at the same time, to adopt a model-independent approach by sharpening those 
theoretical tools that may allow, in the near future, a direct observational test of the effects 
of large-scale magnetic fields on CMB anisotropics. 



54 



References 

[1] Ya. Zeldovich, Sov. Phys. JETP 21, 656 (1965). 

[2] Ya. Zeldovich and I. Novikov, The Structure and Evolution of the Universe, (Chicago 
University Press, Chicago, 1971), Vol.2. 

[3] I. M. Khalatnikov, Sov. Phys. JETP 21, 172 ( 1965). 

[4] E. M. Lifshitz and I. M. Khalatnikov, Sov. Phys. JETP 12, 108 (1960). 

[5] E. M. Lifshitz and I. M. Khalatnikov, Sov. Phys. JETP 12, 558 (1961). 

[6] V. A. Belinskii and I. M. Khalatnikov, Sov. Phys. JETP 30, 1174 (1970). 

[7] M.P. Ryan, and L.C. Shepley, Homogeneous Relativistic Cosmologies, ( Princeton Uni- 
versity Press, Princeton 1975). 

[8] J. D. Barrow, Phys. Rev. D 55, 7451 (1997). 

[9] J. D. Barrow, P. G. Ferreira and J. Silk, Phys. Rev. Lett. 78, 3610 (1997). 

[10] M. Giovannini, Phys. Rev. D 59, 123518 (1999). 

[11] J. D. Barrow and R. Maartens, Phys. Rev. D 59, 043502 (1999). 

[12] M. Giovannini, Phys. Rev. D 62, 067301 (2000). 

[13] K. A. Bronnikov, E. N. Chudaeva and G. N. Shikin, Class. Quant. Grav. 21, 3389 
(2004). 

[14] G. Chen et ai, Astrophys. J. 611, 655 (2004). 

[15] R. Durrer, T. Kahniashvili and A. Yates, Phys. Rev. D 58, 123004 (1998). 

[16] K. Subramanian and J. D. Barrow, Phys. Rev. D 58, 083502 (1998). 

[17] K. Subramanian and J. D. Barrow, Phys. Rev. Lett. 81, 3575 (1998). 

[18] T. R. Seshadri and K. Subramanian, Phys. Rev. D 72, 023004 (2005). 

[19] J.Adams, U.Danielsson, D.Grasso, and H. Rubinstein, Phys. Lett. B 388, 253 (1996). 

[20] D. Grasso and H. R. Rubinstein, Phys. Rept. 348, 163 (2001). 

[21] A. Vlasov, Zh. Eksp. Teor. Fiz. 8, 291 (1938); J. Phys. 9 (Moscow), 25 (1945). 

[22] L. D. Landau, J. Phys. (Moscow) 10, 25 (1945). 

55 



[23] E. M. Lifshitz and L. P. Pitaevskii, Physical Kinetics, (Pergamon Press, Oxford, Eng- 
land, 1980). 

[24] D. Biskamp, Non-linear Magnetohydrodynamics (Cambridge University Press, Cam- 
bridge, 1994). 

[25] T. J. M Boyd and J. J. Sanderson, The physics of plasmas, (Cambridge University 
Press, Cambridge, UK, 2003). 

[26] M. Giovannini, Phys. Rev. D 71, 021301 (2005). 

[27] J. Skilling, Phys. of Fluids 14, 2523 (1971). 

[28] S. M. Carroll, G. B. Field and R. Jackiw, Phys. Rev. D 41, 1231 (1990). 

[29] S. M. Carroll, Phys. Rev. Lett. 81, 3067 (1998). 

[30] M. Giovannini, Phys. Rev. D 61, 063004 (2000). 

[31] M. Giovannini, Phys. Rev. D 61, 063502 (2000). 

[32] B. Feng, H. Li, M. z. Li and X. m. Zhang, Phys. Lett. B 620, 27 (2005). 

[33] . K. Sethi and K. Subramanian, Mon. Not. Roy. Astron. Soc. 356, 778 (2005). 

[34] K. R. S. Balaji, R. H. Brandenberger and D. A. Easson, JCAP 0312, 008 (2003). 

[35] A. Lue, L. M. Wang and M. Kamionkowski, Phys. Rev. Lett. 83, 1506 (1999). 

[36] N. F. Lepora, arXiv:gr-qc/9812077l 

[37] A. Kosowsky and A. Loeb, Astrophys. J 461, 1 (1996). 

[38] M. Giovannini, Phys. Rev. D 56, 3198 (1997). 

[39] D. Harari, J. Hayward and M. Zaldarriaga, Phys. Rev. D 55, 1841 (1997). 

[40] M. Giovannini, Int. J. Mod. Phys. D 14, 363 (2005). 

[41] J. J. Sakurai, Modern Quantum Mechanics, (Addison- Wesley, New York, 1985). 

[42] C.-P. Ma and E. Bertschinger, Astrophys. J. 455, 7 (1995). 

[43] W. Hu and N. Sugiyama, Astrophys. J. 444, 489 (1995). 

[44] H. V. Peris et al, Astrophys. J. Suppl. 148, 161 (2003). 

[45] A. Kogut et al, Astrophys. J. Suppl. 148, 213 (2003). 

56 



[46] C. Scoccola, D. Harari, and S. Mollerach, Phys. Rev. D 70, 063003 (2004). 

[47] L. Campanelli, A. D. Dolgov, M. Giannotti and F. L. Villante, Astrophys. J. 616, 1 
(2004). 

[48] A. D. Dolgov, Lectures given at the International School on Astrophysics 'Daniel Cha- 
longe' ( 8th Paris Cosmology Colloquium 2004 and Topical Meeting WMAP and the 
Early Universe, Paris, France, 9-10 Dec 2004), |arXiv:astro-p h/0503447} 



[49] A. Kosowsky, T. Kahniashvili, G. Lavrelashvili and B. Ratra, Phys. Rev. D 71, 043006 
(2005). 



[50] see, for instance, http://www.rssd. esa.int/index.php?project=PLANCK 



[51] see, for instance, http://oberon.romal.infn.it/boomerang/b2k/| 



[52] S. Masi et al, |arXiv:astro-ph/0507509 



[53] see, for instance, http : / / cfcpwork.uchicago.edu /kicp-pr o j ects / quiet / } 



[54] see, for instance, http://www.skads-eu.org/ 



[55] F. Govoni and L. Feretti, Int. J. Mod. Phys. D 13, 1549 (2004). 

[56] B. M. Gaensler, R. Beck and L. Feretti, New Astron. Rev. 48, 1003 (2004). 

[57] N. A. Krall and A. W. Trivelpiece, Principles of Plasma Physics, (San Francisco Press, 
San Francisco 1986). 

[58] D. T. Son, Phys. Rev. D 59, 063008 (1999). 

[59] A. Brandenburg, K. Enqvist and P. Olesen Phys. Rev. D 54, 1291 (1996). 

[60] A. Brandenburg, K. Enqvist and P. Olesen Phys. Lett. B 392 , 395 (1997). 

[61] P. Olesen, Phys. Lett. B 398, 321 (1997). 

[62] K. Enqvist, Int. J. Mod. Phys. D 7, 331 (1998). 

[63] M. Giovannini, Int. J. Mod. Phys. D 13, 391 (2004). 

[64] L. M. Widrow, Rev. Mod. Phys. 74, 775 (2003). 

[65] L. Pogosian, T. Vachaspati, and S. Winitzki Phys. Rev. D 65, 083502 (2002). 

[66] M. Giovannini and M. E. Shaposhnikov, Phys. Rev. Lett. 80, 22 (1998). 

[67] M. Giovannini and M. E. Shaposhnikov, Phys. Rev. D 57, 2186 (1998). 

57 



[68] T. Kahniashvili and B. Ratra, Phys. Rev. D 71, 103006 (2005). 

[69] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions (Dover, New 
York, 1972). 

[70] I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, Series and Products (fifth 
edition), (Academic Press, New York, 1994). 

[71] M. Giovannini and M. E. Shaposhnikov, Phys. Rev. D 62, 103512 (2000). 

[72] A. Mack, T. Kahniashvili and A. Kosowsky, Phys. Rev. D 65, 123004 (2002). 

[73] K. Kajantie, et al. Phys. Rev. Lett. 77, 2887 (1996). 

[74] T. Kibble and A. Vilenkin, Phys. Rev. D 52, 679 (1995). 

[75] J. Ahonen and K. Enqvist, Phys. Rev. D 57, 664 (1997). 

[76] G. Baym, D. Bodeker and L. McLerran, Phys. Rev. D 53, 662 (1996). 

[77] J. Quashnock, A. Loeb and D. Spergel, Astrophys. J. Lett. 344, L49 (1989). 

[78] B. Cheng and A. Olinto, Phys. Rev. D 50, 2421 (1994). 

[79] G. Sigl, A. Olinto, and K. Jedamzik, Phys. Rev. D55, 4852 (1997). 

[80] A. D. Dolgov, |hep-ph/0110293t 



[81] M. Giovannini, Phys. Rev. D 62, 123505 (2000). 

[82] M. Giovannini, Phys. Rev. D 64, 061301 (2001). 

[83] K. Bamba and J. Yokoyama, Phys. Rev. D 70, 083508 (2004). 

[84] O. Bertolami and R. Monteiro, Phys. Rev. D 71, 123525 (2005). 

[85] K. E. Kunze, Phys. Lett. B 623, 1 (2005). 

[86] M. Gasperini, M. Giovannini and G. Veneziano, Phys. Rev. Lett. 75, 3796 (1995). 

[87] D. Boyanovsky, M. Simionato and H. J. de Vega, Phys. Rev. D 67, 023502 (2003). 

[88] D. Boyanovsky, H. J. de Vega and M. Simionato, Phys. Rev. D 67, 123505 (2003). 

[89] D. Boyanovsky and H. J. de Vega, |arXiv:astro-ph/0502212l 



[90] K. Enqvist, A. Jokinen and A. Mazumdar, JCAP 0411, 001 (2004). 



58 



[91] A. Lazarian, Astron. Astrophys. 264, 326 (1992). 
[92] M. J. Rees, Lect. Notes Phys. 664, 1 (2005). 

[93] A. Lazarian, E. Vishniac and J. Cho, Astrophys. J. 603, 180 (2004); Lect. Notes Phys. 
614, 376 (2003). 

[94] A. Brandenburg and K. Subramanian, Phys. Kept. 417, 1 (2005). 
[95] R. Kulsrud, Annu. Rev. Astron. Astrophys. 37, 37 (1999). 
[96] R. Kulsrud and S. Anderson, Astrophys. J. 396, 606 (1992). 
[97] I. Brown and R. Crittenden, Phys. Rev. D 72, 063002 (2005). 
[98] J. M. Bardeen, Phys. Rev. D 22, 1882 (1980). 

[99] C. G. Tsagas and J. D. Barrow, Class. Quant. Grav. 14, 2539 (1997). 
[100] C. Tsagas and R. Maartens, Phys. Rev. D 61, 083519 (2000). 
[101] C. G. Tsagas, Class. Quant. Grav. 22, 393 (2005) 
[102] M. Giovannini, Phys. Rev. D 58, 124027 (1998). 

[103] T. Kahniashvili, G. Gogoberidze and B. Ratra, |arXi v:astro-ph/0505628| 

[104] C. G. Tsagas, Class. Quant. Grav. 19, 3709 (2002). 

[105] T. J. Battefeld and R. Brandenberger, Phys. Rev. D 70, 121302 (2004). 

[106] M. Giovannini, Phys. Rev. D 70, 103509 (2004). 

[107] M. Giovannini, Class. Quant. Grav. 22, 363 (2005). 

[108] K . Subramanian and J. D. Barrow, Mon. Not. Roy Astron. Soc. 335, L57 (2002). 

[109] K. Subramanian, T. R. Seshadri and J. D. Barrow, Mon. Not. Roy. Astron. Soc. 344, 
L31 (2003). 

[110] A. Lewis, Phys. Rev. D 70, 043011 (2004). 
[Ill] A. Lewis, Phys. Rev. D 70, 043518 (2004). 
[112] M. Giovannini, Phys. Lett. B 622, 349 (2005). 
[113] M. Giovannini, Phys. Rev. D 70, 123507 (2004). 



59 



[114 

[lis; 

[116 
[117; 

[lis; 

[119 
[120 

[121 

[122; 



S. Koh and C. H. Lee, Phys. Rev. D 62,083509 (2000). 

M. Gasperini, M. Giovannini and G. Veneziano, Phys. Rev. D 52, 6651 (1995). 

C. A. Clarkson, A. A. Coley, R. Maartens and C. G. Tsagas, Class. Quant. Grav. 20, 
1519 (2003). 

K. Jedamzik, V. Katalinic and A. V. Olinto, Phys. Rev. Lett. 85, 700 (2000). 



A. Zizzo and C. Burigana, arXiv:astro-ph/0505259| 



P. Dineen and P. Coles, Mon. Not. Roy. Astron. Soc. 347, 52 (2004). 

G. Bernard!, E. Carretti, R. Fabbri, C. Sbarra and S. Cortiglioni, 



|arXiv:astro-ph/0508278[ 



P. D. Naselsky, A. G. Doroshkevich and O. V. Verkhodanov, Astrophys. J. 599, L53 
(2003). 

P. D. Naselsky, L-Y. Chiang, P. Olesen, O. Verkhodanov, Astrophys. J. 615, 45 (2004). 



60 



